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Abstract: A major ingredient in Higgs searches at the Tevatron and LHC is the elimination 
of backgrounds with jets. In current H — )■ WW — )• luiu searches, jet algorithms are used 
to veto central jets to obtain a 0-jet sample, which is then analyzed to discover the Higgs 
signal. Imposing this tight jet veto induces large double logarithms which significantly modify 
the Higgs production cross section. These jet- veto logarithms are presently only accounted 
for at fixed order or with the leading-logarithmic summation from parton-shower Monte 
Carlos. Here we consider Higgs production with an inclusive event-shape variable for the 
jet veto, namely beam thrust T^nu which has a close correspondence with a traditional pT 
jet veto. Tcra allows us to systematically sum the large jet-veto logarithms to higher orders 
and to provide better estimates for theoretical uncertainties. We present results for the 0-jet 
Higgs production cross section from gluon fusion at next-to-next-to-leading-logarithmic order 
(NNLL), fully incorporating fixed-order results at next-to-next-to- leading order (NNLO). At 
this order the scale uncertainty is 15 — 20%, depending on the cut, implying that a larger 
scale uncertainty should be used in current Tevatron bounds on the Higgs. 
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1 Introduction 



The discovery of the Higgs boson is a major goal of the Large Hadron Collider (LHC) and 
current analyses at the Tevatron. The decay H — >• WW* is the dominant channel for Higgs 
masses niH ^ 130 GeV. Hence, the H — >• W~^W~ — )• l'^v^~i' channel has strong discovery 
potential and plays a very important role for early searches that are statistically limited. It 
is the dominant channel in the current Tevatron exclusion limit [1]. The presence of the 
final-state neutrinos does not allow the reconstruction of the Higgs invariant mass, and hence 
sideband methods cannot be used for this channel to determine the backgrounds directly 
from data. At the LHC and Tevatron, ti — )• W^W~bb events constitute a large background, 
dominating the signal by a factor of 10 to 40 depending on the Higgs mass and center-of-mass 
energy. Requiring a minimum missing energy is not effective against this background since 
it also contains two neutrinos. To eliminate the huge background from top-quark decays one 
imposes a stringent jet veto to define a 0-jet sample for the search, where one only allows soft 
jets with < p™*. The latest ATLAS study [2] vetoes any jets with transverse momentum 
Pj,^ > 20 GeV and pseudorapidity {if^^l < 4.8, which reduces the tt background to a negligible 
level. The latest CMS study [3] rejects all events that have jets with pip > 25 GeV and 
{if^^l < 2.5, which reduces this background by a factor of ~ 40. After the jet veto, the main 
irreducible background stems from the direct production channel — ?• M^"*"!^", which at this 
point still dominates the signal by a factor of about 4:1. The final discrimination against 
this and other backgrounds is achieved by exploiting several kinematic variables [4]. 

The Tevatron Higgs searches analyze their data using a jet algorithm and Monte Carlo to 
implement a jet veto and divide the data into 0-jet, 1-jet, and > 2-jet samples for all jets with 
Pt^ ^ 15 GeV and \r]^^^\ < 2.4 — 2.5 [1, 5, 6]. For mn > 130 GeV the sensitivity is completely 
dominated by the 0-jet and 1-jet samples in — )• WW. At lower Higgs masses, the WH, 
ZH, and vector-boson-fusion production channels with higher jet multiplicities are included 
to increase sensitivity. With the latest update from ICHEP 2010 [7], the Tevatron excludes a 
range of Higgs masses mn = 158 — 175 GeV at 95% confidence level. For these exclusion limits 
it is important to have a good theoretical understanding of the jet production cross sections 
and a reliable estimate of theory uncertainties separately for each jet bin, as emphasized in 
ref. [8]. The theory uncertainties in the Higgs production cross section were investigated 
recently in refs. [9-11]. For their 0-jet bin, the Tevatron analyses use an uncertainty of 7%, 
which is taken from the fixed next-to-next-to-leading order (NNLO) analysis of the 0-jet bin 
in ref. [8]. With our resummed next-to-next-to-leading logarithmic order (NNLL) plus NNLO 
calculation of a 0-jet cross section we will see that the perturbative uncertainties are actually 
larger, ~ 20%, due to the presence of large logarithms that are not accounted for in the 
fixed-order analysis. 

Theoretically, the inclusive Higgs production cross section has been studied extensively in 
the literature and is known to NNLO [12-19] and including NLO electroweak corrections [20- 
22] (for reviews and additional references see e.g. refs. [23, 24]). However, Higgs production 
in a 0-jet sample differs substantially from inclusive Higgs production. In particular, the jet 
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veto induces large double logarithms a" ln™(p^*/7n//) with m < 2n that are not present 
in the inclusive cross section, and also induces a sizable dependence on the choice of jet 
algorithm used to define the veto (see e.g. ref. [25]). Theoretical studies of the jet veto 
are available in fixed-order calculations at NNLO [26-28], and include additional kinematic 
selection cuts [8, 25, 29, 30] (see also ref. [31]). 

Currently, the only method available to experiments to incorporate the effect of the jet 
veto and the accompanying large logarithms beyond fixed order is to use parton-shower Monte 
Carlos, such as MC@NLO [32, 33], POWHEG [34, 35], Pythia [36, 37], and Herwig [38, 
39]. This allows one to take into account the dependence of the 0-jet sample on the choice of 
jet algorithm, but for the large logarithms it limits the accuracy to the leading-logarithmic 
summation provided by the parton shower. The comparison [8, 25, 28] of the results at fixed 
NLO with those from MC@NLO, Herwig, and Pythia (the latter two reweighted to the 
total NLO cross section), shows differences of 20 — 30%, cf. tables 4 and 1 of refs. [8, 25] 
respectively. This shows the importance of resumming the phase-space logarithms caused 
by the jet veto. Furthermore, the Herwig and Pythia parton-level results obtained in 
ref. [8] differ by about 15%, which is an indication that subleading phase-space logarithms 
are relevant. 

Theoretically, one can also study the Higgs production as a function of the Higgs trans- 
verse momentum, p^, both in fixed-order perturbation theory for large [40-43] and with 
a resummation of logarithms of p!f at small p^ [44-51]. A further method is the so-called 
joint factorization [48, 52], which allows one to simultaneously resum logarithms at threshold 
and small p^ by introducing pT-dependent PDFs. For H — >• W~^W~ — ?• i~^i'i~iy the missing 
neutrino momenta make a direct measurement of small p^ impossible. Instead the NNLL 
resummed p^f spectrum [53] is used to reweight the Pythia Higgs spectrum in the Tevatron 
search, which is important for estimating the efficiency of selection cuts. The study of pT- 
resummation is also motivated by the fact that the jet veto automatically forces pff to be 
small, see e.g. refs. [8, 25, 54]. However, the logarithms at small summed at NNLL differ 
from those induced by the jet veto. Thus studies of the small-p-^ spectrum can only provide 
a qualitative template for the effect of the jet veto.^ 

In this paper we explore a jet veto in pp — t- HX at the LHC and pp — ?• HX at the Tevatron 
using an inclusive kinematic variable called beam thrust [56]. Beam thrust does not require 
a jet algorithm and is well-suited for carrying out higher-order logarithmic resummation. 
It allows us to directly predict a 0-jet Higgs production cross section using factorization 
techniques without relying on parton showers or hadronization models from Monte Carlo. We 
will present results for both the differential beam-thrust spectrum as well as the integrated 
pp ^ H + Oj cross section with a cut on beam thrust working at NNLL and including the 
NNLO corrections.^ With the large logarithms under control, we are also able to perform a 

^On the other hand, the hadronic Et spectrum could be considered for a central jet veto, and the resum- 
mation at small Et was carried out in ref. [55] at NLL order. 

^Our NNLL resummation is in the jet-veto variable, and is not the same as NNLL threshold resummations 
for the total cross section [57-62]. 
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realistic assessment of the perturbative theory uncertainties. Since a factorization theorem 
exists for the beam thrust spectrum we are also able to rigorously account for the leading 
effect of nonperturbative hadronization corrections. A final advantage of beam thrust is that 
the cross section for the dominant irreducible background, pp — )• VFVF + Oj, can be computed 
with precisely the same jet veto and similar precision, which we leave to future work. 

While H — >• WW provides the most obvious motivation for studying the effect of jet 
vetoes, one can also consider the case oi H ^ 77. Here, the Higgs signal appears as a 
small bump in the 77 invariant mass spectrum on top of a smooth but overwhelming QCD 
background. The signal and background are separated from each other by a combined fit 
to both. The main reducible backgrounds are pp — t- jj and pp — t- j7, while the irreducible 
background comes from QCD diphoton production, pp — t- 77. Experimentally, it is still 
advantageous to separate the data into 0-jet, 1-jet, and > 2-jet samples because in each sample 
the background has a different shape, which helps to gain sensitivity in the fit. However, this 
separation introduces the same theoretical issues as for the jet veto in i7 — )■ WW. Beam 
thrust provides a continuous measure of the 0-jettiness of an event. Hence, instead of using 
separate jet samples it may be useful to perform a combined fit to the beam thrust and 
77 invariant mass spectra. The theoretical formulas presented here can be used to study 
i7 — )• 77, and we briefly comment on this, however we choose to focus on — )• WW. 

In i7 — )• WW, where missing energy plays an important role, the appropriate version of 
beam thrust is defined in the hadronic center-of-mass frame [56, 63] by 

r = ^ , %m = Yl \PkT\ e"l'"=l =^iEk- \pl\) ■ (1.1) 

The central jet veto using beam thrust is implemented by requiring Tcm '^H ^ CQuivalently 
r <C 1. Since the mass of the Higgs, mn, is unknown, for our analysis the dimension-one 
variable Tcm is more convenient than the dimensionless r. The sum over k in eq. (1.1) runs 
over all particles in the final state, excluding the signal leptons from the W decays. Here pkT 
and rjk are the measured transverse momentum and rapidity of particle k with respect to the 
beam axis (taken to be the z axis).'^ For simplicity we assume all particles to be massless. 

To see that a cut on Tcm ^ fnn vetoes central jets, first note that the absolute value 
in eq. (1.1) divides all particles k into two hemispheres r]k,pl. > and Vk^Pk < 0- We can 
now distinguish between energetic particles with Ek ~ mn and soft particles with Ek <^ mn. 
The latter only give small contributions to Tcm- Energetic particles moving in the forward 
direction have Ef^ — \p^\ <^ mn, so they also contribute only small amounts. In particular, 
unmeasured particles beyond the rapidity reach of the detector are exponentially suppressed, 

give negligible contributions. On the other hand, energetic 
particles in the central region have -B^ ~ IPfcl ~ ~ and give a large contribution. 

^Just as for jet algorithms, experimentally the sum will be over pseudo-particles constructed from calorime- 
ter clusters and possibly supplemented by tracking information. Using information from the tracking systems 
is important to reduce the impact of pile-up as it allows to distinguish particles originating from the primary 
hard interaction from those due to secondary minimum-bias interactions. 
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Figure 1. Comparison of different relations between p^"' and 7^™' for the NNLO cross section, 
where the left panel is for the Tevatron and the right panel is for the LHC. The relation 7^™* ~ 
rauip'T^ /raH)^ yields the same leading large logarithm at 0{as) and also the best overall agreement 
at NNLO. Here we used MSTW2008 NNLO PDFs [64] and evaluate the cross section at ^ = ruH. 



Therefore, a cut Tcm < ^ ''^h provides an inclusive veto on central energetic jets 

without requiring a jet algorithm. 

An important question is how the 0-jet cross section cr(7^™*) with the jet-veto cut Tcm. ^ 
7^'^^* compares to the more standard cr(p™*) with a traditional jet-veto cut on the maximum 
PT of the jets, p'^ < p™*. To relate the uncertainties due to the large logarithms in these 
two cross sections, we can compare their leading double- logarithmic terms at 0{as) using our 
computation and the one in ref. [26]: 

a(7;r)oc(l-^hi2ZS^ + ...), a(p5Vt)oc(l-^^ln2^ + ...). (1.2) 

To obtain agreement for the leading-logarithmic terms in eq. (1.2) the correct correspondence 
between the two variables is 



„cut \ 

7S«-™„^^) . (1.3) 

We can check the accuracy of this relation for the two jet vetoes at NNLO numerically 
using the FEHiP program [27, 43] This fixed-order comparison contains not only leading 
logarithms, but also subleading logarithms and non-logarithmic terms. The ratio of the 
NNLO cross sections using different trial relations for the correspondence between 7^™* and 
are shown in figure 1. With the relation in eq. (1.3) the NNLO cross sections differ by 
< 2% at the Tevatron, and < 7% at the LHC, throughout the range of interesting cuts. If we 
multiply the prefactor in eq. (1.3) by a factor of 1/2 or 2 then the agreement is substantially 
worse, close to that of the dotted and dashed curves in figure 1. This confirms that eq. (1.3) 
provides a realistic estimate for the correspondence. However, it does not directly test the 
correspondence for the cross sections with resummation at NLL order or beyond. 
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Figure 2. Comparison of tlie Higgs signal and tt background using Pythia. The differential spectrum 
in %ni is shown on the left, and in p™'^'^, the px of the hardest jet, on the right. For the jet algorithm 
we use the anti-fct algorithm with R = 0.4, only considering jets with |?7j°'| < 2.5 or jr/j'^'l < 4.8. 

To illustrate the relative size of the H — )• WW signal compared to the tt — )• WWbb 
background as a function of either Tcm or the pT of the hardest jet, p^'^^, we use Pythia 
8 [37] to simulate gg ^ H ^ WW for mn = 165 GeV and tt WWbb events. In both 
cases we turn off multiple interactions in Pythia, since the corresponding uncertainty is 
hard to estimate without dedicated LHC tunes. Following the selection cuts from ATLAS 
in ref. [2] we force one W to decay into an electron and one into a muon. We then require 
both leptons to have pT > 15 GeV and |r/| < 2.5. For the dilepton invariant mass we require 
12 GeV < mei < 300 GeV, and for the missing transverse momentum, p^^^^ > 30 GeV. We 
have not attempted to implement any lepton isolation criteria since they should have a similar 
effect on the Higgs signal and tt background. For the pT jet veto we define jets using the 
anti-fcj algorithm [65] with R = 0.4 implemented in the Fast Jet package [66]. The results for 
the differential cross section in and p^^^ after the above cuts are shown in figure 2, where 
the normalization corresponds to the total cross sections (Jgg^H = 8pb and a^l = 163 pb (see 
e.g. ref. [67]). Note that the above selection cuts have no effect on the shape of the Higgs 
signal and a small 5 — 20% effect on the shape of the ti background. In this simulation a 
signal to background ratio of one is achieved with cuts Ton < 31 GeV, p^'"^^ < 32 GeV for 
\r]\ < 2.5, and p™*^^ < 33 GeV for \rj\ < 4.8. It will be very interesting to see the performance 
of Tcm in a full experimental analysis including a 6-jet veto from 6-tagging which will further 
improve the suppression of t — )• Wb decays with only small effects on the Higgs signal. 

We have also tested the correspondence between the Tcm and p^* variables using partonic 
Pythia 8 Higgs samples for the LHC at 7TeV. The cut px < is apphed for R = 
0.4 anti-fcr jets with rapidities |?/| < r/'^"*. For r/'^"* = 4.8 the variable correspondence is 
roughly midway between Tm = p'r^ and the relation in eq. (1.3), whereas for r/'^"* = 2.5 the 
correspondence is closer to Tm = p™*- For the Tevatron the correspondence is also closer 
to Tm = P™^ with less dependence on ry™*. To estimate the impact of our results on the 
uncertainties for the p™* jet veto we will consider the range between eq. (1.3) and Tm = Pj^^- 
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Further discussion on how to apply our results to the experimental analyses using reweighted 
Monte Carlo samples is left to section 4. 

Including the resummation of large logarithms for Tcm ^ fnH, the production cross 
section from gluon fusion, gg — > H, is given by the factorization theorem [56] 



da 

dT^ii, 

^55 



(Jo Hgg{mt,mH,fi) j*^^ j d*a dife Bg{ta,Xa,fi) Bg{th, Xb, fi) 

dcr"'^ 
dTcm 



,Si'(%^-'-^^±^.,)+^. (1.4) 



where 



m 

Xa = e , Xb = e , ctq = ^ (l-^J 

Ecm is the total center-of-mass energy, and Y is the rapidity of the Higgs.^ The limits on the 
Y integration are ln(m/f/£'cm) ^Y<— lii{mH / Ecm)- 

In this paper we focus our attention on the Higgs production cross section. The leptonic 
decay of the Higgs does not alter the factorization structure for the summation of large 
logarithms in the first term in eq. (1.4), where it can be included straightforwardly as was 
done in ref . [56] for the simpler case oi pp ^ Z/^ ^ . Its effect on the second term can 
be more involved. Including the Higgs decay is of course important in practical applications, 
which use additional leptonic variables to discriminate against the pp — )■ WW background. 
A further investigation of these effects using factorization is left for future work. 

By using a cut on Tcm < 7^m' to implement the jet veto, the resulting large double 
logarithms in the 0-jet cross section have the form a" ln'"(7^'^*/m//) with m <2n. Measuring 
Tcra introduces two new energy scales into the problem. In addition to the hard-interaction 
scale ~ rnn^ one is now sensitive to an intermediate beam scale /i^ ~ TcnjnH and 
a soft scale ^5 ~ T^m- In the first term in eq. (1.4), the physics at each of these scales 
is factorized into separate hard, beam, and soft functions, Hgg^ Bg, , which are briefly 
discussed below. The veto induced logarithms are systematically summed using this factorized 
result for the singular terms in the cross section. These functions and the nonsingular cross 
section components, da'^^/dTcm, are discussed in detail in section 2. The full expression in 
eq. (1.4) applies for any value of Tcm, and reduces to the fixed-order result when Tcm — i^H- 

When Tcm ^ ^'t-H the absence of additional hard jets in the final state implies that the 
dominant corrections appearing at fin are hard virtual corrections, which are described by 
the hard function, Hgg{mt,m'jj, hh)- It contains the virtual top-quark loop that generates 
the effective ggH vertex plus the effects of any additionally exchanged hard virtual gluons. 



*For — >■ 77 the Higgs rapidity Y is measurable. With no additional jets in the event it provides the 
boost of the partonic hard collision relative to the hadronic center-of-mass frame. In this case one can account 
for this boost in the definition of beam thrust, Tb = \PkT\ e^'*""^^', which effectively defines beam thrust 
in the partonic center-of-mass frame. Just as for Tcm, a jet veto is obtained by imposing 7s rriH- The 
factorization theorem for the gluon- fusion production cross section for Tb "C mn is the same as in eq. (1.4) 
but with Y set to zero inside S'^ [56]. The difference between da/Tm and dcr/dTs first appears at NLO and 
NNLL and is numerically small, at the 4% level. 
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The jet veto explicitly restricts the energetic initial-state radiation (ISR) emitted by the 
incoming gluon to be colhnear to the proton direction. As a result, the energetic ISR cannot 
be described by the evolution of the standard parton distribution functions (PDFs), which 
would treat it fully inclusively. In this situation, as discussed in detail in refs. [56, 68], the 
initial state containing the colliding gluon is described by a gluon beam function, Bg{t, x,^b), 
which depends on the momentum fraction x and spacelike virtuality — t < of the gluon 
annihilated in the hard interaction. The beam function can be computed as [68, 69] 



Here, fj{^,fiB) is the standard PDF describing the probability to find a parton j with light- 
cone momentum fraction in the proton, which is probed at the beam scale /ib- The virtual 
and real collinear ISR emitted by the parton j builds up an incoming jet and is described by 
the perturbative coefficients Igj{t,x/^, (Ib)- At tree level, a gluon from the proton directly 
enters the hard interaction without emitting any radiation, so Bg{t,x, ^ib) = ^{'t)fg{^it^B)- 
Beyond tree level, real emissions decrease the parton's momentum fraction to a; < ,^ and push 
it off shell with —t < 0. Tcra for small values is given by 

Ton. = ^ + ^ + + OiT^l) , (1.7) 

ruH mn 

where the ta- and t^-dependent terms are the total contributions from forward and backward 
collinear ISR. Here 7^^^* is the total contribution from soft radiation and is determined by 
the beam-thrust soft function S^j^ {T^^\ ^s)- Neither T;^''* nor tafi are physical observables 
that can be measured separately. Hence, in eq. (1.4) we integrate over ta and th subject to the 
constraint in eq. (1.7), where the integration limits are determined by ta^i) > and 7^^^* > 0. 

In section 2 we describe all the ingredients required for our calculation of the 0-jet Higgs 
production cross section from gluon fusion at NNLL-I-NNLO. The hard, beam, and soft func- 
tions are discussed in sections 2.1, 2.2, and 2.3, respectively. In sections 2.4 and 2.5 we describe 
how we add the nonsingular NNLO corrections, which are terms not contained in the NNLL 
result. The treatment of running renormalization scales is described in section 2.6, the impact 
of vr^ summation and PDF choices in section 2.7, and the size of hadronization corrections 
in section 2.8. Details of the calculations are relegated to appendices. (In appendix A we 
calculate the one-loop matching of the gluon beam function onto gluon and quark PDFs, and 
verify at one loop that the IR divergences of the gluon beam function match those of the gluon 
PDF. In appendix B we present analytic fixed-order results for the hard and beam functions 
with terms up to NNLO, as well as results for the singular NLO and NNLO beam thrust 
cross section.) In section 3 we present our results for the Higgs production cross section as a 
function of beam thrust up to NNLL-I-NNLO order. In section 3.1 we study the convergence 
of our resummed predictions. In section 3.2 we compare our resummed to the fixed-order pre- 
dictions, and our main results for the theoretical scale uncertainties are presented in figs. 15 
and 16. The origin of the large X-factors for Higgs production is discussed in section 3.3. 
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Section 4 contains our conclusions and outlook, including comments on the implications of 
our results for the current Tevatron Higgs limits. Readers not interested in technical details 
should focus their reading on the introduction to section 2 (skipping its subsections), and 
then read sections 3 and 4. 



2 Components of the Calculation 

The differential cross section for Tcm in eq. (1.4) can be separated into a singular and nonsin- 
gular piece 

dfj / do-" . dcr'"'\r. . .^/AqcdV 



+ 



l + O 



■niH J 



(2.1) 



dTcm V dTcm dTcm ^ 

Including the renormalization group running of the hard, beam, and soft functions, we have 
aoHgg{mt,m'jj,f_iH)UHimH,IJ.H,tJ-) JdY Jdtadth (2.2) 

dt'aBg{ta - t'a,Xa,fiB) C/^ (t^, /ifi , /i) j dt'f^Bg{tb - t'b,Xb,flB) f/^ (4, , ^) 



dTcr 



[ dk S'^J (Tcm - - — ^° ^'^ - k, us) Usik, iis,lj) 
J V ruH ) 



Equation (2.2) is valid to all orders in perturbation theory and is derived in ref. [56] using 
the formalism of soft-collinear effective theory (SCET) [70-74]. In addition we will consider 
the cumulant, 

o{rZ')= r^^^dTcm;^, (2.3) 

Jo CI /cm 

which gives the cross section with the jet-veto cut Tcxa < T^m*- For cr(T^™*) the relevant scales 
are hh ^ mn, - VmiT^H, and fis - Vi^^. 

Letting ?; — iO be the Fourier conjugate variable to r = Tcm/i^H, the Fourier-transformed 
singular cross section exponentiates and has the form 

In—— ~ Invlaslnv)'' + (aslnv)^ + as(aslnv)^ H , (2.4) 

dv 

where k > 1. The three sets of terms represent the LL, NLL, and NNLL corrections, respec- 
tively. As usual for problems involving Sudakov double logarithms, the summation happens 
in the exponent of the cross section, which sums a much larger set of terms compared to 
counting the leading logarithms in the cross section. To sum the terms in eq. (2.4) to all 
orders in a^, the hard function, Hgg, beam functions, Bg, and soft function, , in eq. (2.2) 
are each evaluated at their natural scales, fin — rnu-, fJ-B — VTcm^H, fJ-s — Tcm, and are then 
evolved to the common scale n by their respective renormalization group evolution factors 
Uh, U^, and Us to sum the series of large logarithms. In table 1 we show various orders 
in resummed perturbation theory and the corresponding accuracy needed for the matching 
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3-loop 


NNLO 


NNLL'+NNLO 


NNLO 


NNLO 


2-loop 


3-loop 


3-loop 


NNLO 


N^LL+NNLO 


NNLO 


NNLO 


3-loop 


4-loop 


4-loop 


NNLO 



Table 1. The order counting we use in fixed-order and rcsummcd perturbation theory. The last two 
rows are beyond the level of our calculations here, but are discussed in the text. 



(i.e. the fixed-order results for the hard, beam, and soft functions) and anomalous dimensions 
(7x, Lcusp) that enter the singular corrections. To NNLL order we require the NLO fixed- 
order corrections for Hgg^ Bg, and , as well as the two-loop non-cusp and three-loop cusp 
anomalous dimensions in the evolution factors, and the three-loop running of q^. 

The nonsingular contributions, da'^^/dTcm in eq. (2.1), are 0{Tcm/mH) suppressed rela- 
tive to the resummed contribution, da^/dTcm- They become important at large Tcm and are 
required to ensure that the resummed results also reproduce the fixed-order cross section at 
a given order. 

For the various combinations in table 1 we show the order at which nonsingular corrections 
are included, which for consistency agrees with the order for the singular matching corrections. 
For example, to include the fixed NLO corrections in the NLL result requires including both 
the singular and nonsingular NLO terms, which we denote as NLL'-I-NLO. Similarly at one 
higher order we would obtain NNLL'-I-NNLO. The prime in both cases refers to the fact that 
the matching corrections in the resummed result are included at one higher order than what 
would be necessary for the resummation only. The complete NNLO matching corrections for 
the beam and soft functions, which we would need at NNLL' and N'^LL, are not available 
at present. Instead, for our final result, which we denote as NNLL-I-NNLO, we only include 
the /Li-dependent NNLO terms in Hgg, Bg, and , which we compute using the two- loop 
RGEs. The remaining /U-independent NNLO terms are added in addition to the nonsingular 
NNLO terms, as discussed in section 2.5, such that the fixed-order expansion of our final 
result always reproduces the complete NNLO expression. 

In the following sections 2.1 to 2.3, the hard, beam, and soft function are discussed in 
turn, including expressions for their fixed-order corrections as well as their NNLL evolution. 
The one-loop results for the hard and soft function are easily obtained from known results. 
The one-loop calculation for the gluon beam function is performed in appendix A. 3. The 
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anomalous dimensions are all known and given in appendix B.3. The basic SCET ingredients 
relevant to our context are reviewed in refs. [56, 68]. To obtain numerical results for the cross 
section, we use the identities from App. B of ref. [75] to evaluate the required convolutions of 
the various plus distributions in the fixed-order results and evolution kernels. In section 2.4 
we discuss how to extract the nonsingular contributions at NLO and NNLO, and in section 2.5 
how these are combined with the resummed singular result to give our final result valid to 
NNLL+NNLO. 

The scale /x in eq. (2.2) is an arbitrary auxiliary scale and the cross section is manifestly 
independent of it at each order in resummed perturbation theory. This fact can be used 
to eliminate one of the evolution factors by setting /x equal to one of /i/^, fiB, or fis- The 
relevant factorization scales in the resummed result at which a fixed-order perturbative series 
is evaluated are the three scales fin, fJ-B, and fis- Hence, their dependence only cancels out 
up to the order one is working, and the residual dependence on these scales can be used to 
provide an improved estimate of theoretical uncertainties from higher orders in perturbation 
theory. The choice of scales used for our central value and to estimate the perturbative 
uncertainties is discussed in section 2.6. Finally, in section 2.7 we briefly discuss the effect 
that the vr^ summation and the order of the used PDFs have on our results. 

2.1 Hard Virtual Corrections 

The hard function contains hard virtual corrections at the scale of order tuh, including the 
virtual top-quark loop that generates the effective ggH vertex. It is obtained by matching 
the full theory onto the effective Hamiltonian in SCET 

^eff = E / d'^idc^s C,,H{mt, 2h ■ 62) (26i • 62)5/.. • (2-5) 

ni,n2 

Here, H denotes the physical Higgs field and v = {\^Gf)~^^'^ = 246 GeV the Higgs vacuum 
expectation value. The Bn,uj fields are gauge invariant fields in SCET that describe energetic 
gluons with large momentum bi = ooini/2, where ni are unit light-cone vectors, = 0. The 
matching coefficient CggH depends on the top-quark mass and the invariant mass 2bi ■ 62 of 
the two gluons. For the case we are interested in we have 2bi • 62 = q'^, where q is the total 
momentum of the Higgs, i.e. of the WW or 77 pair. In addition to the operator shown in 
eq. (2.5), there are also operators where the Higgs couples to two collinear quark fields. The 
tree-level matching onto these operators is proportional to the light quark mass, m^, and are 
numerically very small. There are potentially larger matching contributions from QCD loops 
where the Higgs couples to a top quark, but these are also mq/niH suppressed due to helicity 
conservation. Hence, we neglect these collinear quark operators in our analysis. 
The hard function is defined as 

Hgg{mt,q'^,fM) = \CggH{mt,q'^,^i)\'^ . (2.6) 

It is evaluated at q^ = in eq. (2.2) because we consider the production of an on-shell 
Higgs. (Including the decay of the Higgs, the cross section differential in q'^ is proportional 
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to aoLHgg{mt, q'^ , fi), where L contains the squared Higgs propagator and decay matrix ele- 
ment. In the narrow width approximation L reduces to L = 6{q'^ — m|^)Br, where Br is the 
appropriate Higgs branching ratio, e.g. Br(ff — t- WW) or Br(i7 — t- 77).) 

By matching onto eq. (2.5) we integrate out ah degrees of freedom above the scale fin^ 
which are the heavy top quark as well as gluons and light quarks with offshellness above fin- 
This can be done in either one or two steps. In the one-step matching used here we integrate 
out both the top quark and hard off-shell modes at the same time. This allows us to keep 
the full dependence on m'jj/m^. In pure dimensional regularization with MS the matching 
coefficient Cggnifnt, q'^ , fJ-n) is given by the infrared-finite part of the full mt-dependent ggH 
form factor, which is known analytically at NLO (corresponding to two loops) [76, 77] and in 
an expansion m q'^/m^ at NNLO (three loops) [78, 79]. 

We write the Wilson coefficient as 



V Am 



1 + 



iO 



+ 



«^(^h) 

(4^)2 



iO 



Ami 



+ 



V4m| 



(2.7) 



where F(°)(0) = 1. At NNLL we need the NLO coefficients 



In XH + — 

D 



fW(0) = 5Ca-3Cf- 



(2. 



The dependence of F^'^\z) and F^^\z) on z = q'^/{Amf ), which encodes the mt dependence, 
is given in eq. (B.l). At NNLL+NNLO we also need to include the NNLO terms that depend 
logarithmically on the hard scale fin, which follow from the two- loop RGE of the Wilson 
coefficient (see eq. (B.12)), and are given by 



C^^HxH, z) = hj\ In^ xu + ^Ca/3o In^ xh + Ca 



4 7r2\„ 5 



In^ Xu 



+ 



(2.9) 



The remaining ^//-independent NNLO terms are contained in F^'^^i^z). Although these are 
known in an expansion in z, we do not include them, since the corresponding ^-independent 
NNLO terms are not known for the beam and soft functions. 

To minimize the large logarithms in Cggn we should evaluate eq. (2.7) at the hard scale 
\XH with ^ q^ ^ ml- For the simplest choice /u^ = q^ the double logarithms of —q^/n^ 
are not minimized since they give rise to additional vr^ terms from the analytic continuation 
of the form factor from spacelike to timelike argument, ln^(— 1 — iO) = — vr^, which causes 
rather large perturbative corrections. These vr^ terms can be summed along with the double 
logarithms by taking /_f// = — iy^Q^ or in our case = — i^^H [80-83]. For Higgs production 
this method was applied in refs. [84, 85], where it was shown to improve the perturbative 
convergence of the hard matching coefficient. Starting at NNLO, the expansion of CggH 
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contains single logarithms lii{m^/fi^), which in eq. (2.7) are contained as InxH in C(2) with 
a compensating — ln(— 42; — iO) in F^'^\z), which are not large since mu/mt ~ 1. In eq. (2.7), 
Usi^n) is defined for nj = 5 flavors. When written in terms of asifin) with nj = 6 flavors 
similar ln(m^/mti|^) terms would already appear at NLO. The additional terms that are 
induced by using an imaginary scale in these logarithms are small, because the imaginary 
part of as{—iniH) is much smaller than its real part. 

The alternative two-step matching is briefly discussed in appendix B.l, where we compare 
results with the literature. In this case, one first integrates out the top quark at the scale 
rrit and then matches QCD onto SCET at the slightly lower scale fin — ttih- This allows 
one to sum the logarithms of rriH/mt at the expense of neglecting mj^/m^ corrections. Since 
parametrically mH/mt — 1, we use the one-step matching above. Note that we do not include 
electroweak corrections whose predominant effect (of order 5%) is on the normalization of the 
cross section through the hard function [20-22, 62, 86, 87]. 

Given the hard matching coefficient at the scale ijlh we use its renormalization group 
evolution to obtain it at any other scale 

Hgg{mt, q^,n) = Hgg{mt,q^,i2H) Uniq^ , fJ-H, fJ-) , (2.10) 
where the evolution factor is given by 

7^ — iO\VHi^J■H,^J■) 2 



Uniq ,/i/f,M) 



KH{m,l^) = -'^Kr(t^H,fJ') + K^flj{fJ'H,fJ-) , r]H{f^H,fJ') = r]rif^H,fJ-) , (2.11) 

and the functions K^inHi 1^)-, Vrif^HTfJ-), and K^{fiH,fi) are given in appendix B.3. They 
vanish for fi = fin and therefore Unio^ , HH, IJ-h) = 1, consistent with eq. (2.10). 

2.2 Gluon Beam Function 

The gluon beam function can be computed in an operator product expansion (OPE) in terms 
of standard gluon and quark PDFs (see appendix A.l for more details), 



Bg{t,X,flB)= 22 J y ■^'MBj/j(^,/^B 



j={g,Q,q} 



(2.12) 



In ref. [69] the Igg matching coefficient was computed at one loop in moment space. The Igq 
and Igq coefficients in the sum over j in eq. (2.12) describe the case where a quark or antiquark 
is taken out of the proton, it radiates a gluon which participates in the hard collision, and 
the quark or antiquark then continues into the final state. These mixing contributions start 
at one loop. Our one-loop calculation of Igj for j = {g, q, q}, which are needed for the gluon 
beam function in the 0-jet Higgs cross section at NNLL, is given in some detail in appendix A 
and follows the analogous computation of the quark beam function in ref. [68]. 
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We write the matching coefficients for the gluon beam function as 

z, /.b) = 6{t) 6{l-z) + {t, z,f^B) + {t, z,f,B), 

Igq{t,Z,flB) = I^gq>{t,Z,HB) + -J^^I^gq'{t,Z,flB) ■ (2.13) 

Our calculation in appendix A yields the one-loop coefficients 

X«(t,z,A.B) = 2CAe{z)\^cJ^)6{l -z) + ^cJJ^)Pgg{z) + smj^YHz)] , 



J^B ^f^B f^B ^t^B' 



X« {t,z,fiB) = 2CFe{z)\^Co(-^) Pg, (z) + 5it) X(i''^) {z)\, (2.14) 



where 



= £i(i - ,)^(L^±^ _ _ - .) , 

T^'/\z) = Pgg{z)ln^ + e{l - z)z . (2.15) 

Here Pgg{z) and Pgq{z) are the g ^ gg and q ^ gq splitting functions given in eq. (A. 16), 
and the Cn{x) denote the standard plus distributions, 



Cn{x) 



ix) In^x 



(2.16) 



defined in eq. (A. 44). From eq. (2.14) we see that the proper scale to evaluate eq. (2.12) is 
/i^ ~ t ~ TcmiTT-H- For our final NNLL+NNLO result we also need the /is-dependent terms 
of the two- loop coefficients, contained in Ig'^ and Igl^ ■ They can be computed from the 
two-loop RGE of the Igj (see eq. (B.6)), which follows from the two-loop RGEs of the beam 
function and the PDFs. Our results for these coefficients are given in appendix B.2. 

Our result for Xgg^ is converted to moment space in eq. (A. 40), and except for a vr^ term, 
agrees with the corresponding moment space result given in eq. (68) of ref . [69] . Another com- 
parison can be made by considering the correspondence with the py-dependent gluon beam 
function from ref. [51], which is given in impact parameter space i/t as Bg{tn,x,yT, fJ-) = 
^ggiin,x/^,yT, IJ-) ® /g(C)^)- Taking the yr ^ limit of their bare result should yield agree- 
ment with our bare beam function. In principle the renormalization could change in this 
limit, but their results indicate that this is not the case. Translating their variable t„ into 
our variables, t„ = t/z, the li'niy,^^Qlgg{tn, z,yT, fJ-) from ref. [51] agrees with our result in 
eq. (2.14). In ref. [88] the authors changed their variable definition from t„ to our t = tnZ.^ 
Ref. [88] also calculates the other lij coefficients at one loop. Our result for the mixing 
contribution Igq^ disagrees with the yx ^ limit of their Igq. In particular, their constant 



^The translated result for Igg quoted in ref. [88] has a typo, it is missing a term —S{t)Pgg(z) Inz induced 
by rescaling [1/ fj,'^)Co{tn/ lJ.^)Pgg{z). We thank S. Mantry and F. Petriello for confirming this. 
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Figure 3. The gluon beam function integrated up to imax — 0.1(a; 7 TeV)^. The left plot shows 
xi?g (^niax, a;, Mb)- The right plot shows all results relative to the LO result. The solid lines show the 
LO and NLO results with the perturbative uncertainties shown by the bands. The dashed, dotted, 
and dot-dashed lines show the NLO result without quark contribution, in the large x limit, and the 
small X limit, respectively. See the text for further details. 



term is X^q^\z) = —Pgq{z) ln[(l — z)/z] + 2(1 — z)/z which disagrees with ours in eq. (2.15). 
We have also compared the results of ref. [88] for the quark beam function for ut ^ with 
our earlier results in refs. [56, 68]. The coefficient Iq^^ agrees, but the mixing term ijjl^ also 
disagrees. The limyj,_^o result in ref. [88] is missing a term —6{t)Pqg{z) that is present in 
refs. [56, 68]. 

Given the beam function at the scale fiB from eq. (2.12), we can evaluate it at any other 
scale using its renormalization group evolution [68] 



Bg{t,x,fi) = j dt' Bg{t - t',x,i^iB) UBit',HB,tJ') 



(2.17) 



with the evolution kernel 

pKB-lE VB 

UB{t,HB,IJ') 



VB_ 



(4^) + 6{t) 



r(i + m) 

Kb ip,B ,IJ) = ,fJ')+ K^a^ {fiB , ^) , 



r]B{fJ-B,fJ-) = -2?7^(^B,/i) . 



(2.18) 



The plus distribution £^(x) = [9{x)/x^~^]^ is defined in eq. (A. 44), and the functions 
(//_B , /i) , ?7p(/UB, /Li), and K'y{fiB,tJ') are given in appendix B.3. Note that for n = hb 
we have UB{t, hb, IJ-b) = S{t), which is consistent with eq. (2.17). 

To illustrate our results for the gluon beam function we define its integral over t < tmax. 



-Bg(imax, X, ^Lb) = j '^^ ^gi^^ ^' ^(*max " t) ■ 



(2.19) 



In figure 3 we plot -Bg(tmax5 2;, /j-ij) for a representative fixed value of tmax — 

0.1(3; 7TeV)2. 

(Similar plots for the quark and antiquark beam functions can be found in refs. [56, 68].) 
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The left panel shows xBg{t^i^x,x, fis)- The right panel shows the relative corrections to the 
LO result B^°{t^^^,x, i^ib) = fgix,fiB)- We use MSTW2008 NLO PDFs [64] with their 
as{mz) = 0.12018 and two-loop, five-flavor running for q^. The bands show the perturbative 
uncertainties from varying the matching scale fiB- Since at the scale ^ib there are no large 
logarithms in the beam function, the fiB variation can be used as an indicator of higher- 
order perturbative uncertainties. At LO the only scale variation is that of the PDF and the 
minimum and maximum scale variation are obtained by fiB = {\/^max/2, 2\/imax} with fj,B = 
^/tra^LX the central value. At NLO the maximum variation does not occur at the endpoints of 
the range V^max/2 < fJ-B < 2^*^ 



but rather for approximately ^b 



{O.S-v/tmax, 2.0-v/tmax} 



with the central value at //^ = 1-5-v/^max- The ag corrections to the gluon beam function are 
quite large, between 20% to 40%, which is signiflcantly larger than the ~ 10% corrections to 
the quark beam function. The main reason for this is the larger color factor for gluons than 
quarks. 

The size of the various perturbative contributions to the beam function is illustrated 
in figure 3. The dashed line shows the result obtained from Igg, without adding the mixing 
contribution Igg, using the same central value hb = 1.5-v/imax- The mixing contributions are 
only relevant above x > 0.2, and are suppressed at small x, because of their smaller color 
factor compared to Igg and the dominance of the gluon PDF at small x. This means they 
will be numerically small for a light Higgs. 

The dotted line in figure 3 shows the result in the threshold limit (again for /i^ = 
1.5-v/tmax)) where we drop Igg and in addition only keep the terms in Igg that are singular as 
z — 7- 1, and which are expected to become dominant as x — t- 1 in eq. (2.12), 



l'gg'Ht,Z,f^B) 



5{t)5{l 
2 



^I^B 



5{l 



+ 



-L 



( * 



27r 

Coil -z)+ 5{t) 2£i(l -z)- ^6(1 - z) 



vr 



6 



(2.20) 



The dotted line indeed approaches the dashed line for x > 0.2. However, since in this region 
the mixing contributions become important, the threshold result does not provide a good 
approximation to the full result (solid line) anywhere. 

Finally, the dot-dashed line in figure 3 shows the result only keeping the terms singular 
as 2; — )■ (but including tree level). 



l'g^'>{t,Z,flB) 



l'g^%t,Z,flB) 



6{t)6{l 

«s(^b) 

vr 



B) 



IT 



CAe{z)9il-z) 



1 



Co 



r t \1 



m 



Cf9{z)9{1-z) 



f^B 



fl%J Z 



m- 



(2.21) 



which one might expect to dominate for x — )• 0. Since these only have single-logarithmic fx 
dependence, their central value is obtained for hb = V^max- This contribution indeed grows 
towards smaller x, and makes up more than half of the total contribution at x = 0.01, but 
does not yet dominate. 
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2.3 Soft Function 



The soft function S^{k,ns) appearing in eq. (2.2) is defined by the vacuum matrix element 
of a product of Wilson lines. For k ~ Tcm ^ Aqcd, it can be computed in perturbation 
theory. We write the perturbative soft function as 

S^/^,dk,f^s) = S{k) + ^fit^CAS('J{k,f,s) + ^^CA5g)(fc,/.5) + , (2.22) 



where the one- and two-loop coefficients are 
5£)(fc,/.5) = -— a(— )+^5(fc), 

/X5 \fls^ 12 
5g) {k, ^s) = 8Ca -Cs (-) + /3o -£2 (- 

,7 ^2 



+ 



9 



9 12 



1 



4 8^2 
3 + — 
/ k 



Ca + -/3o 



Af5 



-A 



k 



99 



(2.23) 



In ref. [56] the quark beam-thrust soft function was obtained from the one-loop hemisphere 
soft function for outgoing jets [89, 90]. The gluon beam-thrust soft function has Wilson 
lines in the adjoint rather than fundamental representation and at one loop Sp^j.^ is obtained 
from the quark result by simply replacing Cp hy Ca- The //g-dependent terms needed at 
NNLL-I-NNLO are shown in eq. (2.23) and were obtained by perturbatively solving the two- 
loop RGE of the soft function (see eq. (B.13)). The determination of the /i5-independent 

('2 (5) 

constant term, Sgg S{k), requires the two-loop calculation of the soft function. 

The RG evolution of the soft function has the same structure as that of the beam function. 



5|f (fc, /i) = / dk' S'j^{k - k', f,s) Us{k', t,s,f^) 



(2.24) 



with the evolution kernel 
Us{k,fis,l^) 



r(i + ris) 



+ sik) 



(2.25) 



where £^(x) = [6{x)/x^~'^']+ is defined in eq. (A. 44), and Kj,{fis, fJ.), r?p(/X5,/x), and Kj{^s,fJ') 
are given in appendix B.3. 

The nonperturbative corrections can be modeled and included using the methods of 
refs. [75, 91]. The perturbative component 3^^,.^ and nonperturbative component F^^ of the 
soft function can be factorized as 



S'j!{k,f^s) = Jdk'S'pUk-k',f,s)F^^{k') 



(2.26) 



At small Tcm ~ Aqcd the nonperturbative corrections to the soft function are important. 
When the spectrum is dominated by perturbative momenta with Tcm ^ Aqcd eq. (2.26) can 
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be expanded in an OPE as 



5f = SZdk,,s) - 2^)f + , (2.27) 



dA; V fc3 ; ' 

where the leading power correction is determined by the dimension-one nonperturbative pa- 
rameter Of = / dk' {k' /2)F39{k') which is parametrically ©(Aqcd)- The positivity of F39{k) 
imphes that Qp^ > 0, so the factorization in eq. (2.26) predicts the sign of the correction caused 
by the nonperturbative effects. We will see that this simple OPE result with one nonpertur- 
bative parameter ilf gives an accurate description of the nonperturbative effects in the Tcm 
spectra for the entire region we are interested in, which includes the peak in the distribution. 
The OPE in eq. (2.27) implies that the leading nonperturbative effects can be computed as 
an additive correction to the spectrum 

- P"-^* 2J7f— f^, (2.28) 



and likewise for the cumulant 

^^(T;™* ) = ^ertir^') - 2^^f ^a;ert(7;r ) ■ (2-29) 

To first order in the OPE expansion this is equivalent to a shift in the variable used to eval- 
uate the perturbative spectrum, Tcm %m — 20^, or cumulant, ~^ T^m^ — 2J7f . For 
the cumulant the nonperturbative corrections always reduce the cross section, whereas the 
distribution is reduced before the peak and increased in the tail region. Since the nonsingular 
terms in the cross section are an order of magnitude smaller than the singular terms we can 
also replace by a, that is include the nonsingular da'^^/dTcm in eq. (2.28). For simplicity, 
we will use the purely perturbative result in most of our numerical analysis. However, in sec- 
tion 2.8 we will use eqs. (2.26) and (2.28) to analyze the effect of nonperturbative corrections 
on our predictions. 

2.4 Nonsingular Contributions 

In this section we discuss how we incorporate the nonsingular contributions to the cross section 
using fixed-order perturbation theory. For the two beam thrust cross sections considered in 
this paper, the full cross section in fixed-order perturbation theory can be written as 

da pmfH^L'" ^ 



drdy 



V4m^ 



^ /"^^ VC'ij(^,^,r,y,/i,m/^,mj) /i(^a,At)/i(6,M), (2-30) 

where i,j = g,q,q sum over parton types, and r = Tcm/i^^H-^ To simplify the notation 
in the following we will suppress the rriH and mt dependence of the coefficients Cij. The 



For — !■ 77 where the boost between the partonic and hadronic center-of-mass frames is accounted for 
with tb = Ts/rnH, the appropriate replacements in eq. (2.30) are to take r ^ tb, and in the fourth argument 
0. 
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contributions to the Cjj can be separated into singular and nonsingular parts, 

Cij{Za, Zb, T, Y, /i) = Cij{Za, Zb, T, Y, fl) + CfJ'iZa, Zb, T, Y, fl) , (2.31) 

where the singular terms scale as ~ 1/r modulo logarithms and can be written as 

Ct^iZa, Zb, T, Y, fi) = C-j\za, Zb, Y, ^x) 5{t) + ^ Cfj{Za, Zb, Y, fi) £fc(r) , (2.32) 

A:>0 

where Cki^) = [^{t)0-'^^^^)/'t]+ is defined in eq. (A. 44). The resummed result for the cross 
section in eq. (2.2) sums the singular contributions at small r to all orders, counting (a^ Inr) ~ 
1. The nonsingular contributions, C??, are suppressed relative to the singular ones by 0{t) 
and it suffices to determine them in fixed-order perturbation theory. Hence, we can obtain 
them by simply subtracting the fixed-order expansion of the singular result from the full 
fixed-order result, 

^ = ^-^' ^'-''^ 

and analogously for the cumulant 

^ns,FO(^cut) = CT^O(t^"*) - Ct"'^0(t^"*) . (2.34) 

We must take care to use the same PDFs and renormalization scale fi for both the ex™ and 
^s,FO tgi.j]2s. In our analysis we obtain o-NLO^-^cut-j ^^^^ ^NNLO^^cut^ numerically using the 

publicly available FEHiP program [27, 43], which allows one to obtain the fixed-order NNLO 
Higgs production cross section for generic phase-space cuts. 

At tree level, the only nonzero coefficient is C^^ , and the nonsingular contribution van- 
ishes, £jiis,LO|-^cut-j _ g ]\JLQ^ iT^Q (jk g^j^g nonzero for k < 1 and are fully contained in 
the resummed NNLL result for da^/dr. Hence, we can obtain them by expanding our NNLL 
singular result to fixed next-to-leading order, 

a' {t )=a' (r )|nlo- (.2.35) 

The explicit expressions are given in appendix B.4. Subtracting this from the full NLO result 
we get the nonsingular contribution at NLO, 

^ns,NLO^^cutN ^NLO/^cut\ „s,NNLL/^cut\ | /q QR^ 

cr ' (r ) = cr {t )-a' (r )|nlo- (^-^o) 

In the left panel of figure 4 we plot the NLO nonsingular cross section determined by this 
procedure for three different choices of fi, namely fi = mH/'^,'m}{,'2mH- The scaling of 
the nonsingular distribution implies that it involves only integrable functions, therefore the 
cumulant (j'^^'^^o^^cut^ vanishes for r'^"* — )• 0. The /i dependence of the NLO nonsingular 
cross section is sizeable since this is the leading term in this part of the cross section. This /U 
dependence is canceled by the nonsingular terms at NNLO which we turn to next. 
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Figure 4. The left panel shows the nonsingular contribution to the NLO cross section as a function 
of rs""'™ , for the LHC at 7TeV. The residual NNLO cross section shown in the right panel, is the 
nonsingular NNLO cross section, o-is.nnlo^ pj^g constant term c"^*^^, as explained in the text. 



At NNLO the singular cross section is determined by a result analogous to eq. (2.35) 



s.NNLO/ cut 



.s.NNLL/_cut 



s.NNLOl 



lNNLO,fc>0 

The NNLO singular coefficients C^, are nonzero for — 1 < A; < 3 



(2.37) 



Those for /c > can be 
obtained by expanding the singular NNLL result to fixed NNLO. Their explicit expressions 
are given in appendix B.4. For k = —1 the NNLO contribution to the coefficient C~j^ of 
the 5{t) is not fully contained in the NNLL result. Therefore, the r'^^^-independent k = —1 
contribution is not included in the first term on the right-hand side of eq. (2.37), but in the 
second term. For this second term we proceed as follows. 



First, we write the NNLO contribution to C.,- as 



NNLO 



clj{Za, Zb, Y) + C^AZa, Zb, Y, fj,) + cjfiZa, Zb, Y) 



(2.38) 



(2vr)2 

The first term in brackets denotes the //-independent terms proportional to tt^ that are part 
of the vr^ summation. The second term contains all terms proportional to ln(/i/mj|f), which 
cancel the fi dependence in the NLO result. We can obtain these two contributions analyt- 
ically, and they are given in appendix B.4. The remaining //-independent terms, c]j^, are 
currently not known analytically. They could be obtained when the complete NNLO results 
for the hard, beam, and soft functions become available. Using eq. (2.38), the second term 
on the right-hand side of eq. (2.37) is given by 



a 



s.NNLOl 



c-(/i) + c'^(/i) + c-^(/z). 



(2.39) 



with {x = {7r,/i,res}) 



C^(/i) = fjo 



F(o)( 



m 



H 



V4m^ 



dy 



/ ^ ^ij 



Xa Xb 



(2.40) 
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The /i dependence of c'^(/i) cancels that of o-nS'NLOi-^cut^ terms of O(a^), whereas that 

of c'^(/i) and c^'^^di) only starts at 0{al). Since Xa^b = ("ii//-£'cm)e^'^, the c^(/i) have a 
nontrivial dependence on mn- 

To determine the constant c^^^ numerically, we now consider 



NNLO,A:>0 

Since (jUs.NNLO^^cut^ vanishes as r™* —?■ 0, the coefficient c^"^^ is determined by o-res^^cut-j 
as r'^"* —7- 0, while the nonsingular corrections are given by the remainder (jiff's (^^cut^ _ ^rcs^ 
Hence, we can obtain both by fitting our numerical results for u>^cs^^cut-j different values of 
r'^"* to the following function: 

^rcs^^^ = c'^^s + ao r Inr + ai T + 02 Inr + ag . (2.42) 

The oo through 03 terms are sufficient to describe the r'^*^' dependence of (jii^^'^NLO^^cut-j ^^^^ 
the whole range of r'^"*. The results of the fit for pp collisions at 7TeV and niH = 165 GeV 
for fi = niH, = m/f/2, and fj, = 2m// are shown in the right panel of figure 4. At fj, = mu 
this fit gives 

c''^'^ = 0.86 ± 0.02 , ao = 7.6±0.6, ai = 9.3±1.5, 02 = 3.9 ±1.1, 03 = -9.9 ±1.6. 

(2.43) 

Similarly, for pp collisions at 1.96 TeV, rriH = 165 GeV, and = uih, we obtain 

c'''^" = 0.028 ± 0.001 , ao = 0.28 ±0.03, ai = 0.44 ± 0.08 , 02 = 0.18 ± 0.06 , 

03 = -0.42 ±0.08. (2.44) 

We have checked that we obtain the same values for c^^^ within the uncertainties when the fit 
range is restricted to the region of small r™* . (In this case the fit is not sensitive to 02 and 03 , 
so either one or both of them must be set to zero.) Note that (jNNLO^^cut^ ^s,nnlL|^^^^ 
both diverge as r'^"* — ?• 0. The fact that the difference in (j^cs^^cut^ does not diverge for 
^cut _^ Q pj-Qvides an important cross check between our analytic results for the NNLO 
singular terms and the full numerical NNLO result. The numerical uncertainty from fitting 
c^^^{fi) is much smaller than its ^ dependence, and so can be ignored for our final error analysis. 
The fi dependence of (f°^{fi) comes from the PDFs and the overall (/x), where the latter is by 
far the dominant effect. To see this we can rescale c'^^^(m///2) = 1.35 and d^^{2mH) = 0.58 as 
c^'^^{mH /2)ag{mH) /atirnn /2) = 0.91 and c''®'^(2mi/)af (mi/)/a^(2m//) = 0.83, giving values 
that are very close to the central value d^^{mH) = 0.86. 

Having determined the nonsingular contributions to the cross section we can compare 
their size to the dominant singular terms. In figure 5 we plot the singular, nonsingular, and 
full cross sections at NNLO for fi = mu- The left panel shows the absolute value of these 
components of the differential cross sections (obtained by taking the derivative of (t(t'^"*) with 
respect to r™*). For Tcra ^ t^h the nonsingular terms are an order of magnitude smaller than 



- 21 - 




Figure 5. Comparison of tlic singular, nonsingular, and full cross sections at NNLO for /i = mn- 
The left panel shows the magnitude of the differential cross sections on a logarithmic scale. The right 
panel shows the corresponding cumulant cross sections. 



the singular ones. On the other hand for Tcm ^ itt-h/'^ the singular and nonsingular terms 
become equally important and there is a large cancellation between the two contributions. 
These features of the fixed-order cross section will have implications on our choice of running 
scales discussed in section 2.6. 

To determine the singular NNLO contributions in eq. (2.37) for the above analysis we 
only considered the A: > terms contained in o-^.nnll^ q£ gQ^;^i.gg ^s.nnll ^^^^ contains some 
k = —1 terms at NNLO, in particular the c'^(/i) and c^(/x) contributions, but also parts of 
the contribution from cross terms between the NLO matching corrections. Since we 

know c'^'^(/i) numerically, we are able to determine the missing k = —1 contribution at NNLO 
numerically, which corresponds to the sum of the unknown /u-independent NNLO matching 
corrections to the hard, beam, and soft functions. It is given by the difference 

^s,NNLO ^s,NNLL| i i „rcsr,\ ^s,NNLL| /r. Ar\ 

c{fi) = a' -a' \j^^^Q = c [n) + c^{fi) + c [fJ,) - a ' |NNLO,fc=-i- ^^-^^^ 

Since we include the //-dependent NNLO matching corrections in u'^'^'^^^, its NNLO expan- 
sion is obtained by setting fJ-s = fJ-B = I^H = Thus, we can easily evaluate eq. (2.45) 
numerically. For niH = 165 GeV, we find for the LHC at 7TeV, 

c\'mH/2) = 0.002 , c^niH) = -0.035 , c^lmn) = -0.028 , (2.46) 

and for the Tevatron, 

c''(mH/2) = -0.0043, c\mH) = -0.0026, c^lmn) = -0.0027 . (2.47) 

Comparing this to c^'^^{mH) = 0.86 (LHC) and (f'^^{mH) = 0.028 (Tevatron), we see that these 
coefficients are almost fully accounted for by cross terms between the NLO hard, beam, and 
soft functions. The remaining NNLO terms in are in fact very small, and our NNLL-I-NNLO 
results are therefore numerically very close to the complete NNLL'-I-NNLO result. 
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2.5 Cross Section at NNLL+NNLO 



Using the results of sections 2.1 to 2.4 our final result at NNLL+NNLO for the distribution 
and cumulant is obtained as 



j^NNLL+NNLO 

dTcm 
NNLL+NNLO /^-cut^ 



^^^s,NNLL ^^5 ^^ns,NNLO+7r 
+ + - 



\ ' cr 



a 



dTcm dTcm dTcm 

,s,NNLL/--rcut\ , _<5/T-cut\ , _ns,NNLO+7r2 



( /cm j + ( /cm j + 



(7;r ) . 



(2.48) 



The first term in each equation contains the resummed singular result obtained from eq. (2.2) 
to NNLL order, including the /i-dependent NNLO matching corrections. The last term con- 
tains the NNLO nonsingular corrections determined in the previous subsection, but including 
vr^ summation by using 



da 



ns,NNLO+7r2 



a 



dTcm 

ns, NNLO+TT^/T-cutN 
V 'cm / 



da 



ns,NNLO 



as(/"ns)C. 



dTcm 

ns,NNLO/T-cut\ 
V ' cm / 



vr 



da 



ns,NLOn 



27r 

«s(/^ns) CA 

27r 



dT^i, 



^2^ns,NLO 



(2.49) 



Here UH{m'jj, — i^ns, /^ns) = exp[as(/Xns)Cyl7r/2 + . . .] contains the vr^ summation. There are 
two reasons we include the vr^ summation for the nonsingular terms. First, using SCET 
one can derive factorization theorems for the nonsingular terms when Tcm ^ it^h, and the 
results will involve a combination of leading and subleading hard, jet, and soft functions. 
Many of these terms will have the same LL evolution for their hard functions, and hence they 
predominantly require the same tt^ summation. As a second reason we observe from figure 5 
that there are important cancellations between the singular and nonsingular cross sections for 
T^m ^ f^n/'^- Since the vr^ summation modifies the cross section for all Tcm and T^™* values 
it is important to include it also in the nonsingular terms to not spoil these cancellations. 

The middle terms in eq. (2.48) incorporate the singular NNLO terms that are not repro- 
duced by our resummed NNLL result. At fixed order, they are given by 



dXr 



■cut\ 
cm J 



NNLO 



c^(/x) 



(2.50) 



NNLO 



As we saw in section 2.4, c^{fj.) turns out to be very small numerically, which means we might 
as well neglect the term entirely. We include it for completeness in eq. (2.48), in order to 
formally reproduce the complete NNLO cross section. In fact, at this level other contributions 
that we neglect here, such as the bottom-quark contributions or electroweak corrections, are 
likely more relevant numerically. 

Formally, c^(//) is reproduced by the complete two- loop matching required at NNLL' or 
N^LL. When it is properly incorporated into the predictions at that order it is multiplied by 
a Sudakov exponent that ensures that the total <7(T^j^*) — >■ as T^^* — )• 0. Hence, we can 
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include it in the resummed result by multiplying it with the NNLL evolution factors, 

= C^(/ins) UH{'m'jj,HH,fJ') I dtadthU^{ta,fJ,B,fJ') U^{tb, f^B , fJ-) Us(Tcui " -^^-^j/Us,/! 
d/cm J ^ 

(2.51) 

The scale where we evaluate c^{fJ.) here is an N^LL effect, and so beyond the order we are 
working. We choose fi = fi^g, which is the scale at which we evaluate it'^'^(7^'^"*'* 



' cm 



2.6 Choice of Running Scales 

The factorization theorem in eq. (2.2) resums the singular cross section by evaluating the hard, 
beam, and soft function at their natural scales fin — — fJ,B — Vt'i^h, l^s — titih where 
they have no large logarithms in fixed-order perturbation theory. Their renormalization group 
evolution is then used to connect these functions at a common scale. This resums logarithms 
in the ratios of fin, fJ-B, and /is', which are logarithms of r. The r-spectrum has three distinct 
kinematic regions where this resummation must be handled differently and we will do so using 
r-dependent scales given by profile functions usi^) and UBiT)- Profile functions of this type 
have been previously used to analyze the B — )• spectrum [75] and the thrust event shape 
in e+e~ — )• jets [92]. 

For Aqcb/it^h ^ t ^ 1 the scales hh, IJ'B, and Aqcd are all widely separated 
and the situation is as described above. We define this region to be ri < t < T2. In the 
T < Ti region the scale jis drops below 1 GeV, we have Aqcd/w-h ~ r, and nonperturbative 
corrections to the soft function become important. In this case the scales are ^j^h — —^niH, 
fJ-B - y^^QCDmn, and fis ^ Aqcd- 

Finally for t > T2 we have r ~ 1, the resummation is not important, and the nonsingular 
corrections, which are evaluated at a fixed scale ^Uns) become just as important as the singular 
corrections. In this region there is only one scale ifin = I^B = IJ'B = l^ns — ^H- Furthermore it 
is known from B and thrust [75, 92] that there can be important cancellations between 

the singular and nonsingular terms in this limit, and that to ensure these cancellations take 
place the scales /^^(t) and /^5(t) must converge to \hh\ = /^ns in a region, rather than at 
a single point. To ensure this we make the approach to quadratic for T2 < t < and 
set UBiT) = fJ-si^) = ^IJ-H for T > r3 (recall that ifin > 0). As we saw in section 2.4, the 
singular and nonsingular contributions in our case become equally important for r > 1/2. 
Accordingly, we choose the profile functions such that the scales converge around this value 
and stay equal for larger r. 

A transition between these three regions is given by the following running scales 



;Ub(t) = 1 + CB ^(ts - t)(1 ) \//U ^run(T, fJ-) , 



r 



I-Is{t) = l + es6{T3-T)[l Atmn(T,/U), 

A*ns = A* • (2.52) 
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Figure 6. Profiles for tlic running scales ^h, ^J■B, and i^is- The central lines for fiB and fis show 
our central scale choices. The upper and lower curves for fis and fis correspond to their respective 
variations b) and c) in eq. (2.55). 



For the profile fJ^runiT, ^) we use a combination of two quadratic functions and a linear function 
as in ref. [92]. For t > our choice for /irun(''") A*) ensures that our cross section formula 
becomes precisely the fixed-order result. 



2aT + b 



T2) 



T < Ti , 
Tl < T < T2 
T2 < T < r3 
T > T3, 



/Wo - 



(2.53) 

TI-T2-T3 n-T2-T3 

The expressions for a and b follow from demanding that fJ-rnn{T) is continuous and has a 
continuous derivative. The value of /xq determines the scales at r = 0, while ti^2,3 determine 
the transition between the regions discussed above. For our central value we use the following 
choice of parameters 

5 GeV 



fj. = ruH , eB = es = 0, ^0 = 2 GeV , n 



T2 = 0.4 , T3 = 0.6 . (2.54) 



The corresponding running scales are shown in figure 6. 

Since the factorization theorem is not affected by 0(1) changes of the renormalization 
scales, we should vary them to determine the perturbative uncertainty. For a reasonable 
variation of the above parameters, the cross section is most sensitive to /x, 63 and eg. We 
therefore estimate our uncertainties from higher order terms in perturbation theory by taking 
the envelope of the following three separate variations, 



a) 




= 2^^mH , 


eB 


b) 




= mn , 


eB 


c) 




= ruH , 


eB 



0, 

±0.5, 
0, 



es 
es 
es 



0, 
0, 

±0.5. 



(2.55) 
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The effect of variations b) and c) are shown in figure 6 by the upper and lower curves for fis 
and fj^Sj respectively. The effect of variation a) is to change the overall vertical scale of the 
\i^h\, fJ-B, and fis curves in figure 6 by a factor of 1/2 or 2 as indicated by the arrows. In 
predictions based on fixed-order perturbation theory only a scale variation analogous to a) 
can be considered. 

2.7 PDFs and vr^ Summation 

In this subsection we briefly discuss the choice of the order of the PDFs for our resummed 
results and the effect of the vr^ summation. 

As shown in table 1, by default we use the PDFs that correspond to the order of the match- 
ing corrections, namely LO PDFs at NLL and NLO PDFs at NNLL. Since the MSTW2008 
PDFs [93] are extracted simultaneously with the value of as{mz), by using PDFs at different 
orders we are forced to also use different values of as{mz)- Our NNLL+NNLO results con- 
tain two-loop corrections and so at this order our default is to use the MSTW2008 PDFs [93] 
at NNLO with their as{mz) = 0.11707 and with three- loop, five-flavor running for as{fJ,). 
For our NLL'-|-NLO and NNLL results, which include one-loop matching, we use the cor- 
responding NLO PDFs with their as{mz) = 0.12018 and two-loop, five-flavor running for 
as{iJi). At LL and NLL, which only includes tree-level matching, we use the LO PDFs with 
o-siiT^z) = 0.13939 and one-loop, five-flavor running for as{fi). 

Note that at NLL (NNLL) there is a slight mismatch in the required running of as{fJ.). 
The resummation at this order requires two- loop (three- loop) Q!s(/u) running, whereas the 
used LO (NLO) PDFs employ one- loop (two- loop) running of as(/x). In this case, we use 
the following compromise. We use the appropriate as{mz) and running consistent with the 
PDF set to obtain the numerical value of ag at some required scale. At the same time, in the 
NLL (NNLL) RGE solutions we use the QCD /3 function at the appropriate one higher loop 
order to be consistent with the RGE. There is no such mismatch in the Os running at LL, 
NLL'-I-NLO, and NNLL-I-NNLO, and hence no mismatch for our highest order predictions. 

Various results which test the effect of tt^ resummation and the treatment of PDFs are 
shown for the cumulant cross section, cr(7^™*), for mn = 165 GeV in figure 7. The left panel 
shows the Tevatron case and right panel shows the LHC with Ecm = 7TeV. The lower 
four blue curves show the NLL and the upper two orange curves the NNLL results. (The 
nonsingular corrections do not affect this discussion much, so for simplicity we do not include 
them in this figure.) The solid lines correspond to our default results, while the dashed, 
dot-dashed, and dotted are variations with other choices for the PDFs or vr^ summation. 

As discussed in section 2.1, the hard function contains large a^vr^*" terms, with m < n, 
which in our default results are summed by evaluating the hard function at fin = —iruH- The 
TT^ summation is switched off by taking fiH = mu instead, which is shown by the dotted lines 
in figure 7. The effect of vr^ summation is very large. It almost doubles the NLL cross section 
and increases the NNLL cross section by about 30%. From the fact that the NLL results 
with vr^ summation (blue solid line) is very close to the NNLL result without vr^ summation 
(orange dotted line) we can conclude that the large corrections from NLL to NNLL are caused 
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Figure 7. The effect of tt^ summation and using different orders for the PDFs on the cumulant 
beam thrust cross section for tuh = 165 GeV at the Tevatron (left panel) and the LHC with 7 TeV 
(right panel). Shown are the NLL result with/without tt^ summation and with LO/NLO PDFs, as 
well as the NNLL cross section with/without tt^ summation and NLO PDFs. See the text for further 
explanations. 



by the large vr^ terms in the virtual hard-matching contributions. This result for the cross 
section with a cut on T^m agrees with the observations made in ref. [84] for the total cross 
section. Similarly, we have checked that the vr^ summation in the NNLL result brings it much 
closer to the total NNLO result. As a result, the convergence of the perturbative series is 
significantly improved by including the tt^ summation, and we will make use of this for our 
main results. 

We can also explore the effect of using the NLO PDFs already at NLL, which amounts to 
including some higher-order terms in the NLL result but allows us to use the same value for 
OLskf^z) at NLL and NNLL. The corresponding results are shown by the blue dashed (with 
TT^ summation) and dot-dashed (without vr^ summation) in figure 7. For our default (solid 
blue curve) we use the LO PDFs in the NLL cross section, which increases it by about 10% 
compared to using NLO PDFs at NLL. Thus it moves the NLL result in the direction of the 
NNLL result. The main reason for the upward shift is the higher value of associated with 
the lower-order PDFs, since the Higgs cross section has an overall a^. 

2.8 Nonperturbative Corrections 

As discussed in section 2.3, for Tcm ^ fnu the leading nonperturbative hadronization cor- 
rections to the beam thrust spectrum are given by a nonperturbative soft function F^^ . For 
Tcra ^ Aqcd the dominant nonperturbative effect can be described by an OPE that yields 
a single nonperturbative parameter Vl^^ ~ Aqcdj leading to a shift in the beam thrust spec- 
trum by Tern — ^ 7^m — 2rjf^, and in the cumulant by 7^^* — > T^^^ — 2Qf^. To illustrate the 
size of this nonperturbative effect we consider two values for Vl?^ . First, r^f^ = 0.35 GeV is 
motivated by the fit result for an analogous parameter for dijet quark production in e^e~ col- 
lisions [92]. Second, = 1.0 GeV is motivated by a potential enhancement by Ca/Cf = 9/4 
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Figure 8. Shift to tlie NNLL+NNLO pcrturbative cross section, shown by soUd curves with ilf^ = 0, 
caused by the leading nonperturbative hadronization corrections, shown by the dashed and dotted 
curves for ftf^ = 0.35 GcV and ftf^ = 1.0 GeV, respectively. 

from Casimir scaling for adjoint Wilson lines. This choice also reproduces roughly the size of 
hadronization effects for Higgs production in Pythia. Using eqs. (2.28) and (2.29), results 
from the OPE for the LHC with £"011 = 7TeV are shown in figure 8. Comparing the OPE 
results for the distribution, shown in the left panel, to the full convolution with a model soft 
function using eq. (2.26), we find that the OPE works well for the entire displayed spectrum 
when = 0.35 GeV and for T^m > 10 GeV when Of = 1.0 GeV. (Thus, for Of = 0.35 GeV 
the peak is perturbative.) Examining the right panel of figure 8, we see that at 7^^* — ^0 GeV 
a power correction of Of = 0.35 GeV reduces cr(7^™*) by 3%, while for Of = 1.0 GeV the 
reduction is by 7%. (The results for the Tevatron are very similar, giving reductions by 2% 
and 6%, respectively, for 7^^* = 20 GeV.) The sign of this nonperturbative shift is predicted 
by the factorization theorem, while its magnitude is determined by Of. Examining the Tcm 
spectra from Pythia before and after hadronization, we find that the hadronization cor- 
rection in Pythia is consistent with the nonperturbative shift discussed here with a value 
Of = 1.0 GeV for both the Tevatron and LHC. 

3 Numerical Results 

In this section we present our numerical results for the Higgs production cross section for both 
the differential beam thrust spectrum, da/dTcm, and the cumulant, cr(7^^*), which gives the 
integrated cross section with a cut on beam thrust, Tcm < 'Tcm - We are mostly interested in 
the region of small Tcm or TSa^ which corresponds to the 0-jet region. We will show resummed 
results up to NNLL+NNLO order and also compare with the results obtained in fixed-order 
perturbation theory at NNLO using FEHiP [27, 43]. An explanation of the various orders is 
given at the beginning of section 2 and in table 1. Since our focus in this section is on the 
perturbative results and their uncertainties, we will not include the nonperturbative hadronic 
correction discussed in section 2.8 (i.e. we take Of = 0). 
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The perturbative uncertainties in the resummed predictions are estimated as explained 
in detail in section 2.6. For the fixed-order results we use /i = mH/2 as the default choice 
for the central value, which tends to give a better convergence for the total cross section, 
mimicking the effect of the vr^ summation. The perturbative scale uncertainties at fixed order 
are then evaluated using ^ = vriH and fi = m///4. (We follow ref. [8] and do not vary 
the renormalization and factorization scales independently.) Since our focus here is on the 
perturbative uncertainties, we do not add PDF and as{mz) uncertainties in our plots. We 
have checked that they are essentially independent of the cut on beam thrust and the same 
as for the total inclusive cross section. 

We show results for both the Tevatron and the LHC. For the LHC we always use E^m. = 
7TeV. The results for higher E^ra are qualitatively similar, except for the overall increased 
cross section. For most of our plots we use rriH = 165 GeV, which is near the WW threshold 
and where the current Tevatron limits are most sensitive. We also show some plots that 
illustrate the dependence of our results on mn- 

3.1 Convergence of Resummed Predictions 

To study the convergence of the resummed results, we consider results at three different 
orders: NLL, NLL'+NLO, and NNLL+NNLO, which contain the matching and nonsingular 
corrections at LO, NLO, and NNLO, respectively. We choose NLL instead of LL as our lowest 
order to compare to, since NLL is the lowest order where we get a useful approximation with 
appropriately large scale uncertainties. (The LL results are lower than the NLL ones and also 
have a smaller scale uncertainty, which means that they do not contain enough information 
yet to provide a reasonable lowest-order approximation.) 

In figures 9 and 10 we show the convergence for the differential spectrum and cumulant, 
respectively for the Tevatron (left panels) and the LHC (right panels). In figure 11 we show 
the cumulant for fixed 7^^* as a function of the Higgs mass. We see that the perturbative 
corrections are rather large, as is typical for Higgs production. The convergence within our 
perturbative uncertainty bands is reasonable for the differential spectrum and quite good 
for the cumulant, both for different 7^^* and different mu- The large step from NLL to 
NLL'-I-NLO is mostly due to the NLO matching corrections. As we saw in figure 5, for 
'^cm* ^ "^H the nonsingular terms are much smaller than the singular corrections that we have 
computed analytically. One can also see this by comparing the size of the NLO nonsingular 
terms in the left panel of figure 4 with the full cross section in the right panel of figure 10. 

The beam thrust spectrum in figure 9 is peaked in the 0-jet region at small beam thrust 
Tcra — 5 GeV with a large tail towards higher values. The peak in the spectrum is a perturba- 
tive prediction; hadronization effects only have a mild effect on the peak structure as shown 
above in section 2.8. For the beam thrust spectrum of Drell-Yan, which is Fig. 3 of ref. [94], 
the peak occurs at smaller values, around Tcm ^ 2 GeV, and the tail of the spectrum falls off 
much faster. The reason for the shifted peak and higher tail for Higgs production compared 
to Drell-Yan is that the incoming gluons emit much more initial-state radiation than quarks. 
This is also the main reason why the perturbative uncertainties are still rather large even 
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Figure 9. The beam thrust spectrum for Higgs production for mn = 165 GeV at the Tevatron 
(left) and the LHC for Ecm — 7TeV (right). The bands show the perturbative scale uncertainties as 
explained in section 2.6. 
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Figure 10. Higgs production cross section as a function of 7^™' for rriH = 165 GeV at the Tevatron 
(left) and the LHC with Ecm = 7TeV (right). The bands show the perturbative scale uncertainties as 
explained in section 2.6. 
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Figure 11. Higgs production cross section with a cut on beam thrust as function of mn at the 
Tevatron for T;™* = 10 GeV (left) and the LHC with E„n = 7TeV and T;™* = 20 GeV (right). The 
bands show the perturbative scale uncertainties as explained in section 2.6. 
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Figure 12. Same as figure 10 but without the tt^ summation. 



at the highest order, NNLL+NNLO. One can also see that at the LHC the tail is somewhat 
higher and the peak less pronounced than at the Tevatron. Correspondingly, the cumulant at 
the Tevatron starts to level out earlier than at the LHC. The reason is that due to the higher 
center-of-mass energy at the LHC, more phase space is available for initial-state radiation. 

In figure 12 we illustrate what happens if we turn off the vr^ summation. Comparing 
with figure 10, we see again that the vr^ summation significantly improves the convergence, 
and by reducing the overall size of the fixed-order corrections it also reduces the size of the 
perturbative uncertainties. 

3.2 Comparison of Resummed and Fixed-Order Predictions 

In this subsection, we compare our best resummed result at NNLL+NNLO to the NNLO 
fixed-order prediction without any resummation. In figure 13 we compare both predictions 
for the differential beam-thrust spectrum for the Tevatron (left panel) and the LHC (right 
panel). For small Tcm the large logarithms of Tcral'^^H dominate the cross section. The NNLO 
cross section contains terms up to ln^(7^m/"if/)/'7^m and diverges as 7^m — ^ 0, so we do not 
expect it to provide a good description of the spectrum at small Tcm- In the NNLL+NNLO 
calculation, the series of logarithms is summed to all orders, which regulates the divergences 
and yields a reliable prediction for the cross section. The resummation also enhances the 
radiative tail in the spectrum, because it essentially sums up the effects of multiple emissions 
from ISR. 

In figure 14 we illustrate that the NNLL+NNLO result correctly reproduces the NNLO 
result for large 7^™*. To see this we need to switch off the vr^ summation, because the 
NNLL+NNLO result would otherwise contain higher order vr^ terms at large that are 
absent at fixed NNLO. Furthermore, we use /x = mu for the NNLO central value, and 
/i = 2mii and = muj^ for the NNLO scale uncertainties. In this way, the NNLL+NNLO 
and the NNLO are evaluated at the same scales for large Tcm ^ 0.6 m/f, where the logarith- 
mic resummation is switched off and our running scales satisfy [see section 2.6] p^s^cm) — 
/i^CT^m*) = /"ns = [i-H = "'^H- In figure 14 we see that with these choices the NNLL+NNLO 
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Figure 13. Comparison of tlie beam tlirust spectrum at NNLL+NNLO to tlie fixed NNLO result at 
tlie Tevatron (left) and the LHC with Ecm = 7TeV (right). The bands show the pcrturbative scale 
uncertainties. 
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Figure 14. Illustration that the NNLL+NNLO resummcd result reproduces the fixed NNLO result 
at large beam thrust for the Tevatron (left) and the LHC (right). The resummed result has the tt^ 
summation switched off and /i = mn is used for the central value of the fixed-order result. The bands 
show the pcrturbative uncertainties. See text for further explanations. 

indeed smoothly merges into the NNLO result, including the scale uncertainties, at large 
7^™*, as it should. Examining figure 14 for smaller 7^™* values we see that the resummed 
result starts to deviate from the fixed-order one for 7^"^* < 40 GeV at the Tevatron and 
Vm ^ 50 GeV at the LHC. 

In figures 15 and 16 we compare the full NNLL+NNLO including vr^ summation to 
the NNLO (using again the default fi = ruH/'i, as the central value) for the Tevatron and 
LHC, respectively. The left panels show the cumulant cross section, and the right panels 
show the same results as the relative difference in percent to the central NNLL+NNLO 
curve, which makes it easy to read off uncertainties. The relative plots are cut off below 
"Tcm* — 5 GeV because the resummed cross section goes to zero there. The central value of 
the NNLL+NNLO leaves the fixed-order uncertainty band at 7^™* — 25 GeV at the Tevatron 
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Figure 15. Comparison of tlic NNLL+NNLO result for tire liiggs production cross section as a 
function of 7^™* to the fixed NNLO result for the Tevatron. The bands show the perturbative scale 
uncertainties. The left plot shows the cumulant cross section. The right plot shows the same informa- 
tion as percent difference relative to the NNLL+NNLO central value. 




Figure 16. Same as figure 15 but for the LHC with Ecm = 7TeV. 




Figure 17. Same as the left panels of figures 15 and 16 but plotted up to 7^™' = itlh- 
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Figure 18. Contribution to tlie relative uncertainties in the NNLL+NNLO results shown in figures 15 
and 16 from the individual scale variations. Here, the /i/f, fiB, and fis variations correspond to cases 
a), b), and c) in eq. (2.55) and are shown in figure 6. 



and at 7^^* — 35 GeV at the LHC. Hence, for any lower values the resummation should 
be taken into account. At 7^^* = 20 GeV the central values of the NNLL+NNLO and the 
NNLO already differ by 20% at the Tevatron and over 25% at the LHC, which both quickly 
grows beyond 50% towards 7^*^* = 10 GeV.'^ This clearly shows that it is important to resum 
the higher-order logarithms that are missing in the fixed-order prediction in order to obtain 
reliable predictions in the 0-jet region. This also means that one cannot expect the scale 
variation in the fixed-order result to give a realistic estimate of the size of the missing higher- 
order terms, and hence it should not be used to estimate the perturbative scale uncertainty 
either. In contrast, since the resummation takes into account the presence of large logarithms 
for small 7^™*, we are able to obtain reliable estimates of the perturbative uncertainties. The 
perturbative uncertainties at small 7^™* are larger than those in the fixed-order result, namely 
15 - 20% at NNLL+NNLO for 7;^^^^* = 15-20 GeV. Implications of this for the Higgs search 
are taken up in section 4. 

In figure 17 we show the same comparison of NNLL+NNLO and NNLO cumulants, 
but plotted up to 7^™* = mn- We can see that the central value of the resummed result 
including the vr^ resummation almost exactly reproduces the NNLO result which uses fi = 
mn/'i'- However, the vr^ summation included in the NNLL+NNLO results leads to a reduction 
of the scale uncertainties in the inclusive cross section compared to those at NNLO. For 
the uncertainties at the LHC we find +3% and —5% and for the Tevatron +5% and —9%. 
In figure 18 we show the relative uncertainties at NNLL+NNLO from the individual scale 
variations in eq. (2.55). We can see that at small 7^™* the uncertainties are dominated by 

and ;U5, i.e. variations b) and c) in eq. (2.55). By construction, those variations go to 
zero at large 7^^') where the uncertainties are now completely determined by variation a) in 

^One miglit expect tliat a better agreement could be acliieved by using a dynamical 7^'J,"'-dependent scale 
in the fixed-order prediction. We have checked that using the intermediate scale = in the fixed-order 
result however does not improve its behavior relative to the resummed result, but in fact makes it a bit worse. 
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Figure 19. Comparison of tlic NNLL+NNLO result to tlic fixed NNLO result for the Higgs production 
cross section with a cut on beam thrust as function of mn at the Tevatron for 7^™* = 10 GeV (left) 
and the LHC with Ecm = 7TcV and 7^™* = 20 GeV (right). The bands show the perturbative scale 
uncertainties. 

eq. (2.55) and denoted hh in the figure, which is equivalent to the fixed-order scale variation. 

Finally, figure 19 shows a comparison of NNLL+NNLO to NNLO as a function of the 
Higgs mass for a fixed value of 7^™*. For smaller Higgs masses, the logarithms get smaller 
and the cross section increases, which reduces the relative differences. This does not change 
however our overall conclusions for the importance of the resummation for both the central 
value and determining the perturbative uncertainties. 

3.3 Discussion of ET-Factors 

Using our results we can also address the origin of the large NLO and NNLO K-factors that 
are typically observed for Higgs production. It is sometimes argued that the origin of these 
large X-factors for the inclusive cross section are large perturbative corrections due to hard 
emissions. This is based on the observation that by vetoing hard jets, the fixed-order K- 
factors are reduced. As a result the fixed-order perturbative series in the presence of a jet 
veto actually appears to be better behaved than for the inclusive cross section. 

Our results show that once the large jet-veto logarithms are properly summed, the K- 
factor is mostly independent of the jet veto. Hence, it is not caused by hard real emissions, 
but rather mostly by hard virtual corrections and to a lesser extent by collinear and soft 
virtual and real radiation. In our analysis this can be examined directly by comparing the 
convergence of perturbation theory for the hard, jet, and soft functions. As we have seen in 
section 2.7, and was already observed in ref. [84], by summing the large tt^ terms in the hard 
virtual corrections, the i^-factors are substantially reduced. 

There is a simple reason why the large iC-factors in fixed-order perturbation theory are 
reduced: At NLO, the jet veto reduces the cross section just because it cuts out available 
phase space, and since at LO the cross section is not yet affected by the jet veto, the NLO K- 
factor is reduced accordingly. Essentially, the large negative phase-space logarithms resulting 
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from the jet veto happen to cancel the large positive corrections from hard virtual corrections. 
A similar effect appears at NNLO where the jet veto reduces the available phase space for the 
second jet. Since the hard virtual corrections are independent of the jet veto one can always 
choose a particular value for the jet- veto cut such that they are exactly canceled by the large 
phase-space logarithms at a given order in perturbation theory. However, to conclude that the 
jet veto in general renders the fixed-order perturbative series better behaved one would have 
to know that the same level of cancellation will happen at each order in perturbation theory. 
Since these two types of corrections are a priori not related there is no reason to believe this 
will be the case. Instead, to obtain reliable predictions both types of large corrections should 
be rendered as convergent as possible. For the hard virtual corrections the vr^ summation 
improves the convergence, and for the large phase-space logarithms this is achieved by the 
resummation carried out here. With the resulting cross section we can then obtain more 
realistic estimates for higher-order theoretical uncertainties as discussed in section 3.2. 

4 Conclusions 

A major effort at the LHC and Tevatron is devoted to the search for the Higgs boson. In 
the current Tevatron exclusion limits and the early LHC searches the H — ?■ WW — t- l'^vl~i> 
channel plays an important role, since it is the dominant decay channel for Higgs masses 
above ~ 130 GeV. A large background for this channel are ti — )• WWbb events, which must 
be eliminated by a veto on central jets, so the resulting measurement is pp ^ H + jets. 
Such a jet veto causes large double logarithms in the cross section, which need to be summed 
to all orders in perturbation theory in order to obtain reliable theoretical predictions and 
uncertainties for the H + jet production cross section. 

In this paper we have studied Higgs production from gluon fusion, using the inclusive 
event shape beam thrust, Tcm, to implement a central jet veto. As beam thrust character- 
izes the event as a whole and does not require a jet algorithm, it is well suited to analytic 
calculations. This allows us to resum the jet-veto logarithms to NNLL order, based on the 
factorization theorem for the beam thrust cross section derived in ref. [56]. In our analysis 
we also include the full set of NNLO corrections, such that our final result at NNLL-I-NNLO 
provides the best possible theoretical prediction at any value of beam thrust. 

Our main results are presented in figs. 15 and 16. We find that in the 0-jet region at 
small beam thrust, the resummation of jet-veto logarithms is crucial, and our central value 
at NNLL-I-NNLO for the cross section with a cut on beam thrust, Tcm < T^m^ differs signifi- 
cantly from the fixed-order prediction at NNLO. We also find substantially larger perturbative 
scale uncertainties arising from the jet veto compared to those at NNLO. Since fixed-order 
perturbation theory is not reliable in the presence of large jet-veto logarithms, one also can- 
not expect its scale variation to yield a reliable estimate of perturbative uncertainties due to 
neglecting higher-order corrections. 

At present, the jet-veto logarithms are taken into account in the experimental analyses 
using the leading-logarithmic resummation provided by parton-shower Monte Carlo programs, 
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usually supplemented with some reweighting procedure to reproduce the total NNLO cross 
section. While this might yield a reasonable central value that takes into account the dominant 
effect of the logarithmic resummation, the uncertainties in the so predicted 0-jet cross section 
cannot be taken as those of the inclusive NNLO cross section. They could at best be equal 
to those in our NNLL+NNLO results. In fact, they are probably larger than that, since we 
include the resummation at two orders higher than the LL resummation of parton showers. 
As we have seen in section 3.1, already one order lower, at NLL'+NLO, the perturbative 
uncertainties are much larger than those in our NNLL+NNLO results. 

The conventional method to veto central jets is to use a jet algorithm and require p'^ < 
Pj^^ for all jets in the event. As we saw in section 1, p™* can be related to 7^^* by associating 
Tc^^/niH — {Pj^^ /mn)^ or 7^"^* = p^^, where the former works well at NNLO while the 
latter is favored by Pythia 8. Hence, we can use the perturbative uncertainties in our 
results based on 7^^* as a benchmark for the perturbative uncertainties from large logarithms 
relevant for p^*, on which the current experimental analyses are based. For example, the 
perturbative uncertainties for typical values p™* ~ 20 — 30 GeV at the LHC can be as large as 
those for 7;^^* ~ 10 - 15 GeV, which are 15 - 20% at NNLL+NNLO (for m// = 165 GeV and 
Ecm = 7TeV). In the current Tevatron analyses, the perturbative scale uncertainty in the 
0-jet cross section for = 15 GeV is taken as 7% from the fixed-order analysis in ref. [8]. 
In contrast, the perturbative uncertainties at NNLL+NNLO for values 7^™* < 15 GeV are 
significantly larger, e.g. about 20% for T^^^ = 10 GeV and nin = 165 GeV. In light of 
this, the current Tevatron exclusion limits should be reconsidered with an increased theory 
uncertainty for the 0-jet Higgs cross section. We note that our conclusions about theoretical 
uncertainties are based on a systematic study of the jet veto, and are therefore independent of 
the arguments in ref. [9] proposing the use of a larger range for the fixed-order scale variation. 

To implement our NNLL+NNLO results in the experimental searches, one could reweight 
the partonic beam-thrust spectrum obtained from Monte Carlo to our results. This would 
incorporate both the higher-order resummation of large logarithms in the 0-jet region and 
at the same time the full NNLO result for the total cross section. This also allows a consis- 
tent treatment of the perturbative uncertainties in both regions. To illustrate this, consider 
dividing the total inclusive cross section, atotai, into a 0-jet bin, ctq = (^{Tcm^), and a re- 
maining (> l)-jet bin, a>i = cJtotai — cr(7^m*)- This separation into exclusive jet bins causes 
large logarithms of 7^^"* in both do and a>i which cancel in their sum. This implies that 
the uncertainties due to 7^^* in both bins are anti-correlated. In particular, if one wants to 
consider the theory uncertainties in (Tq ^ind C7>i simultaneously in a Gaussian fashion, one 
has to consider the full theory covariance matrix 



where Aq and A>i are the theory uncertainties of ctq and (t>i, and /9o,>i is their correlation 
coefficient. In the second step above we used the approximation that the uncertainties in (Tq 
and (Ttotal cire uncorrelated and are added in quadrature in a>i. From our results in figure 18 




(4.1) 



-37- 



we can see that this is reasonable since the uncertainties in ao are dominated by the lower 
scales fiB and ns, while those in atotai are determined by hh (which in this case is combined 
in quadrature with the others rather then by taking the envelope). The uncertainty squared 
for the total cross section, do + cr>i, is given by the sum of all entries in the matrix C in 
eq. (4.1). Due to the anti-correlation the uncertainties for the individual jet bins can be larger 
than that for the total cross section. Our numerical results for the theory uncertainties for 
the 0-jet bin directly give Aq. The full correlation can be taken into account by reweighting 
the Monte Carlo to both the central value curve as well as the results obtained from the 
individual scale variations. 

It would also be useful to have a benchmark theoretical uncertainty that is desired for the 
experimental searches, since with further effort our NNLL+NNLO results can be extended 
to NNLL'+NNLO or N^LL+NNLO, which has the potential to reduce the uncertainty in the 
resummed perturbation theory. Given that the theoretical predictions and their uncertainties 
are very sensitive to the jet veto, it would also be useful to implement the jet veto in the ex- 
perimental analyses directly in terms of beam thrust, for which resummed theory predictions 
are available. In addition, a benchmark experimental study of the beam-thrust spectrum can 
be made with Drell-Yan pairs, as advocated in ref. [94]. 

In this paper, we have restricted ourselves to studying the case oi gg ^ H + jets. 
The same methods can be used to calculate the dominant irreducible background from direct 
WW production, i.e., the process pp — )■ WW + jets using beam thrust for the jet veto. 
The generalization of beam thrust to processes with N signal jets is provided by the event 
shape A'^-jettiness introduced in ref. [63]. It can be used in an analogous fashion to study the 
exclusive H + 1 jet and H + 2 jet cross sections, the latter being relevant for Higgs production 
from vector-boson fusion. 
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A NLO Calculation of the Gluon Beam Function 

In this section we compute the NLO gluon beam function. We first recall its definition, 
renormalization, and operator product expansion in terms of standard PDFs and Wilson 
coefficients Xjj, as discussed in detail in ref. [68]. For a discussion of the basic SCET ingredients 
relevant to our context we refer to refs. [56, 68]. We then give the details of the one-loop 
matching calculation, which provides a check on the one-loop anomalous dimension of the 
gluon beam function and yields the Xij at NLO. 

A.l Definition and General Results 

As usual, we use lightlike vectors = and = (1,— n), satisfying "n? 

n ■ n = 2, to write a vector in light-cone coordinates 

= — + p~ — + with p'^ = n ■ p , p~ = n • p . 

The bare gluon beam function operator is defined as 

= -u e{u) Bl^^{Q) 6{ub+ - COP+) [6{u: - Vn)B':^l{0)] ■ (A.2) 

The SCET gluon field describes n-collinear gluons and includes zero-bin subtractions [95]. 
It contains collinear Wilson lines Wn to render it invariant under collinear gauge transforma- 
tions. It is the field after the BPS field redefinition [73] to decouple soft gluons and is thus 
invariant under soft gauge transformations. Hence, 0^^^° is fully gauge invariant. The gluon 
beam function is defined by the proton matrix element of the corresponding renormalized 
operator Og{t,uj, fi), 

Bgit,x = oj/p-,fi) = {pniP^)\Ogit,u;,fi)\pn{P~)) . (A.3) 

The proton states \pn{P~)) have lightlike momentum P^ = P~n^/2, i.e. n is chosen in the 
direction of the proton, and the matrix elements are always implicitly averaged over the 
proton spin. 

In eq. (A.2), p~^ is the momentum operator of the residual plus momentum and acts 
on everything to its right. The additional phase in the position-space operator is included 
to allow us to write the b'^ dependence in terms of 6{u}b~^ — uip'^), as in the second line 
of eq. (A.2). Hence, b~^ measures the total plus momentum of real initial-state radiation, 
i.e. of any intermediate state inserted between the fields. The label operator 5{uj — Vn) in 
eq. (A.2) only acts inside the square brackets and forces the total sum of the minus labels 
of all fields in Bn± to be equal to ui. In eq. (A.3) it sets the fraction of the proton's light- 
cone momentum carried into the hard collision to j; = lo/P~ > 0. At the time of the 
collision, the annihilated collinear gluon propagates in a jet of initial-state radiation, and by 
momentum conservation the small plus momentum of the gluon is —b'^ < 0. Hence, it is 



= = 0, 
(A.l) 
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spacelike and —t = w(— 6+) < measures its transverse virtuality. The operator Og{t,Lj) is 
RPI-III invariant, because the transformation of the overall u compensates that of 6{u} — Vn)- 
Therefore, the gluon beam function only depends on the RPI-III invariant variables x and t. 

The RGE of the beam function follows from the renormalization of Og{t,oj, fi). To all 
orders in perturbation theory, it takes the form 

d f 

11 ^Bg{t, X,fl)= J dt' 7|(i - t', /i) Bg{t', X, fi) , (A.4) 

with the anomalous dimension 

7|,(i,/i) = -2r?,3p[a,(^)] ^Co(^) +7|,[«s(/u)]5(t) . (A.5) 



cusp L"-*s Vr*/J 2 \ 2 

Here, Co{x) = [6{x)/x]^ is defined in eq. (A. 44), rcusp(as) is the gluon cusp anomalous 
dimension, and 7^(0!^) denotes the non-cusp part, whose coefficients up to three loops are 
given in appendix B. The RGE of the beam function is identical to that of the jet function. 
It sums double logarithms of t/fj? and unlike the PDF evolution does not change the x 
dependence and does not mix different parton types. The solution of eq. (A.4) was already 
given in eq. (2.17). 

The standard PDFs are defined in SCET by the matrix elements [74] 

/i(e = Uj'/p-,fl) = {pn{P^)\Qi{uj',fl)\pn{P^)) , (A.6) 

where Qi{u}' , fx) are the MS renormalized operators corresponding to the bare operators 
Q^-^(a;0 = e(a;0xn(O)|[<5(^' -P„)xn(0)] , 

Qbare(^/) ^ 0(^') tr{ |xn(0) [5(a.' - P„)xn(0)] } . (A.7) 

The 9{u}') here separates the quark and antiquark PDFs and is included to keep analogous 
definitions for the PDFs and beam functions. 

Compared to eq. (A.7), the fields in the beam- function operator in eq. (A. 2) are separated 
along the n light-cone with large separation y~ corresponding to the small momentum 6+ <^ uj. 
Hence, by performing an operator product expansion about the limit y~ — )• we can expand 
the beam-function operator Og{t,io, fi) in terms of a sum over Qj(w', /u), 

Og{t,co,n) = Y^ ['^lJt,^,^Q^{J,f,) + o(y-), (A.8) 

— J UJ \ id y \ OJ ^ 

j 

where the functional form of the matching coefficients Iij{t, z, fj,) is again determined by RPI- 
III invariance. The proton matrix element of eq. (A.8) yields the OPE for the gluon beam 
function quoted in eq. (1.6), 

Bg{t, X,fl)=^l^ Igj (t, |, ^) (e, /i) [1 + O (^) 



(A.9) 
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The power corrections scale like (AQ(-,j-)/a;)/fe+, where t = a;6+ and Aq^^q/o; is the typical plus 
momentum of the partons in the proton, and involve higher-twist proton structure functions. 
For t ~ A^cD the OPE would require an infinite set of higher-twist proton structure functions, 
which means the beam functions essentially become nonperturbative ^'''-dependent PDFs. On 
the other hand, for t S> ^qcd^ ^^'^ compute the matching coefficients lij in eq. (A. 9) in 
perturbation theory at the beam scale fj,^ ~ t. 

The matching calculation is carried out as usual by taking partonic quark and gluon 
matrix elements of the operators on both sides of eq. (A. 8), which we denote as 

Bg/j{t,Z = Uj/p^,l^l) = {jn{P^)\Og{t,U}, fi)\jn{p^)) , 

fi/j(.z = w/p^,/^) = (jn(P~)|Qi(w,^)|i„(p^)) , (A. 10) 

and analogously for the bare quantities. Here, \jn{p~)) is an n-collinear gluon or quark state, 
\gn{p~)) or \qnip~)), with momentum = p^n^ /2 and p^ > 0, and we average the matrix 
elements over both color and spin. We denote the partonic momentum fractions hy z = uj/p~ 
to distinguish them from the hadronic x and The operators are normalized such that at 
tree level we simply have 

Bf/^it, z, ^) = 5gj 5{t) 5{l - z) , ^) = 6,, 6{l - z) , (A.ll) 

which yields 

llf{t,z,fi) = 6gj6{t)6{l-z), Bl^\t,x,fi)=6{t)fg{x,fi). (A.12) 

Beyond tree level x and ^ will be different, because the momentum fraction ^ provided by the 
PDFs is modified by the collinear radiation encoded in the Igj{t, z, fi). 

Expanding eq. (A. 8) to NLO and using the above tree-level results we find 

B^^lit, z, fi) = X« (t, z, I,) + 6it)fll (z, fi) , 

=^lHt,z,^^)+6{t)fl]l{z,f^). (A.13) 

Hence, to compute the one-loop coefficients 2^^^(i, z, //) we have to compute the one-loop 
gluon and quark matrix elements of the gluon operators Og{t, uj, fi) and Qg{uj, fj,). (The quark 
and antiquark coefficients are identical, Igq = Igq, so we do not need to consider separate 
matrix elements with external antiquarks.) The corresponding one-loop calculation for the 
quark operators to determine the matching coefficients Iqj for the quark beam function is 
given in detail in ref. [68]. There, the calculation is performed in two different ways, first using 
different regulators for IR and UV, and second using dimensional regularization for both IR 
and UV. In the first case, the UV and IR divergences can be separated allowing one to obtain 
the UV renormalization of the operators and at the same time check that the IR divergences in 
the beam function match those in the PDFs, so the matching coefficients are IR finite. In the 
second case, the calculation is much quicker, because most diagrams are scaleless and vanish. 
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but it does not allow one to distinguish the IR and UV divergences. This means one can only 
check that the IR divergences cancel in the matching if the renormalization of both beam 
function and PDF are already known. Alternatively, one can assume that the IR divergences 
cancel and use the known renormalization of the PDFs to obtain the renormalization of the 
beam function. In our case, the renormalization of the gluon beam function is already known 
from the general analysis of ref. [68], so we only perform the second, simpler, calculation, 
regulating both UV and IR in dimensional regularization using the MS scheme. 



A.2 The Gluon PDF at One Loop 

As explained in ref. [68], the definition in eq. (A. 6) is equivalent to the standard definition 
of the PDFs in QCD. In particular, the collinear fields in eq. (A. 7) do not require zero-bin 
subtractions, because the soft region does not contribute to the PDFs. Therefore, we can use 
the standard MS renormalization of the PDFs, 



(A.14) 



where j = {q,q,g}, and the entries in the matrix z(-[z,ii) are a series of 1/e poles with 
coefficients in terms of the renormalized MS coupling as{^). At one loop with d = 4 — 2e, we 
have for the standard gluon PDF [96], 



/.) = 8{l -z) + e{z) [CaW + ^/?o 5(1 - z) 

zUz,f^) = -^^^^^eiz)p,,iz), 



(A.15) 



where /3o = {IICa — 4^^/T/)/3, is the lowest order coefficient of the QCD /3 function, and the 
g ^ gg and q ^ gq splitting functions are 



Pggiz) = 2Coil-z)z + 29{l-z) 



rl 



+ z{l-z) 



= 2e{i-z) 

Pgg{z)=e{l-Z) 



a-z) + 

l + (l-z)2 



1 ~ -2 n' 

+ + z{l-z) 



(A.16) 



with Co{l - z) = [9{1 - z)/{l - z)]+ defined in eq. (A.44). 
Expanding eq. (A.14) to one loop, we get 



bare (1) , 



{z) = zj;';'\z,^,)+r'{z,^.). 



(A.17) 



In pure dimensional regularization, the only Lorentz invariant quantity fgj'j^iz) can depend 
on is z. Since z is dimensionless, all loop diagrams contributing to the PDF calculation 
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y-/uj,-b+ 



y-/uj,-b+ 



p-£,a \ 




(d) 




(e) 




(f) 



Figure 20. One- loop Feynman diagrams for the gluon beam function. The minus momentum w is 
incoming at the vertex and the momentum is outgoing. Graphs (b), (d), and (e) have symmetric 
counterparts, which are equal to the ones shown. The corresponding diagrams with the external lines 
crossed or involving a four-gluon vertex vanish and are not shown. 



vanish, f^'J^j^^^\z) = 0, meaning that the IR and UV divergences in f^'J^j^^^\z) are equal with 
opposite signs. Thus, from eqs. (A. 17) and (A. 15) we get 



].?^e{z)[CAPgg{z) + \PoK^-z) 



1 as{fi)CF 



e 2-K 

where the e poles are IR divergences. 



0{z)P,g{z)^ 



(A.18) 



A. 3 The Gluon Beam Function at One Loop 

We now turn to the one-loop calculation of the gluon beam function. The relevant one-loop 
diagrams are shown in figure 20, with the Feynman rules given below in eq. (A. 21). Regulating 
both the UV and IR in dimensional regularization, the virtual diagram in figure 20(d) and 
its symmetric counterpart vanish because there is only the p~ momentum of the external 
gluon flowing into the loop, which is insufficient to give a nonzero Lorentz-invariant result 
for the loop integral. The diagram in figure 20(e), corresponding to the contribution from 
wave-function renormalization, vanishes for the same reason. Thus, in the on-shell scheme, 
both the wave-function renormalization constant as well as the residue i?^ entering the 
LSZ formula are equal to one. (A different scheme would give contributions to both and 
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p,a,cx q,b,f3 



p,a,<y. 

(a) (b) 
Figure 21. The SCET Feynman rules for the gluon field strength at 0{g°) (a) and 0{g) (b). 



that cancel each other in the final result.) The diagram in figure 20(c) vanishes, because 
the external gluons have perpendicular polarization, so contracting the two in eq. (A. 21) 
leads to n • n = in the numerator. 

Hence, the only diagrams we have to compute are those in figures 20(a) and 20(b) (plus 
its symmetric counterpart), which determine B^J^^°^^\t, z), and figure 20(f), which deter- 
mines B^J^^^^^\t, z). Since the diagrams have real radiation in the intermediate state, we can 
compute them as the discontinuity of the matrix element of the corresponding time-ordered 
product of fields [68] 



sj/7(t,^/p-) 



(tj)Discj>o 
X ^j„(p-)|r{B^ 



^ giij/-/(2a;) 



_ n 



[6iLo-Vn)B'::[iO)]}\Up-))- (A.19) 



For our calculation we need the Feynman rules for the SCET gluon field strength 



with one and two external gluons. They are illustrated in figure 21 and are given by 



figure 21(a) =5^'^(<7r 
figure 21(b) = igr' 



a -a 

p'j_n 
n-p 



(0) ^la 



ip) 



(A.20) 



(A.21) 



+ 



fp±_ 






\n-q 


n-p) 


n-{p + q)_ 



^9r'B^r!l'''^{p,q)^ 



n-p n-q 

where p and q are both taken as incoming. We abbreviate the triple gluon vertex as 

gr'^W'^iPi - P2y + 9'''{P2 - PzT + g'^ips - PiY] = 9r''vr''{pi,p2,P3) , (A.22) 

where all momenta are incoming and momentum conservation holds, P1+P2+P3 = 0. Finally, 
figure 20(e) requires the collinear quark-gluon vertex, which we write as 



n-p n-q n-pn-q J I I 



(A.23) 
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where p and q are the momenta of the outgoing and incoming quark hues, respectively, so 
the gluon carries incoming momentum p — q. 

Let us start with the diagram in figure 20(a). The average over the polarizations of the 
external gluons with momentum p^ = gives 

J—ye-e*^ = -4^, (A.24) 

pol 

where d = 4 — 2e. We also average over the color of the external gluons. This choice for the 
polarizations and colors of the external states is just a matter of calculational convenience. 
Using states with fixed polarization and color gives identical results. To make our expressions 
more palatable, we first perform the color algebra, which yields 

jaecfbde ^fc ^fd ^ _ (^25) 



iV2-l 

The diagram is then given by 
figure 20(a) 

./e^-/.2y52c.^ ^ /■ d'^-¥x (5 + 4/z2)^l + (l + lA)t 

= .(.)^(1 - .)2 [1^ + .(1 - .) + |] r(l + (^) ^ . 

In the second line we integrated and t~ and introduced a Feynman parameter a. The 
last line follows from the standard integrals and the discontinuities listed in appendix A. 4. 
Expanding in e, using eq. (A. 45), we arrive at 



figure 20(a) = "^^^ ^ g(z)g(l-z)2 
27r 



-+z 1-z +- 



z ' ' 2 

(A.27) 

The color algebra for the diagram in figure 20(b) yields 

___^aec ^/c^/e6 ^ (^.28) 
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Including a factor of two for its mirror graph, we get 



figure 20(b) = -2i 



5^(7^(9(0;) Disc 



d'^i 5{i--uj)6{i++b+) 
{27ry (£2 + i0)[(£-p)2 + i0] 



X 1^ (p, e - p, -I) Bf^ ^^{i) B^^l\^(p - A -p) 



d-2 
as{n)CA 



e{z) 



l + z 



(2vr)2 

e{z){i + z)z 



1-z 



Discj>o / da 



d 



1 



using again the relations listed in appendix A. 4. Expanding in e yields 



hgure 20(b) = 0[z) 



--6it) + ^Co(^ 



—5{l-z) + Co{l-z)z{l + z) 



12- 



TT 



+ ^£1 [ — )6il-z) + 6{t)z{l + z) £1(1 - z) - £0(1 - 2) In z - — (^(1 



12 



(A.30) 



Note that crossing the external lines of these diagrams corresponds to p — )• —p (the 
polarization and color are symmetric and hence unaffected). By changing i ^ —I this yields 
the same as the original diagram with the replacement (5(£~— a;)(5(£"'"+6"'") — )• b{t~^ijj)b{l'^—h^') 
which leads to 

(A.31) 



a 



A= 1 -- 1 + 



a 



in the denominator of the loop integrals. Since a, z > these graphs have no discontinuity 
for t > and vanish. Other one-loop diagrams that could in principle connect the two gluon 
fields, e.g. involving a four gluon vertex, also vanish, because they require to cut one of the 
external lines entering the gluon operator, which sets A = t in the loop integral and leads to 
a vanishing discontinuity for t > 0. 

Adding up the two nonvanishing diagrams in figures 20(a) and 20(b), we find 



^bare(l)^^ 



z, = ^iM^e(z) 



27r 
2 



2 1 



r£o 



b(\-z)-U{t)P,,{z) 



t 



1 



+ -o A 5(1 - ^) + 



+ 5(t) £1(1 -z) 



2{\-z + z') 



2\2 



TT 



P,g{z)\xiz-—b(\-z) 



(A.32) 



where Pgg{z) is the g ^ gg splitting function given in eq. (A. 16). 

The mixing contribution of the quark PDF into the gluon beam function originates from 
the diagram in figure 20(f). The spin and color averages are 



(A.33) 



spins 
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The diagram is thus given by 
figure 20(f) 

^■7K ,,2. f 

^^Cf ^(w)Disc; 



47r 



i>0 



d'^£ {p-- e-)5{i--uj)5{£++b+) 



(2vr) 



An J (27r)2 
27r 



+ -:7 Disc 



-z)r(l + e)(e^^/i2)^ 



't>o 



'2+i0)2[(p-^)2 + i0] 
) 



(A.34) 



da (1 — a 



€Z 



sinvre 6'(t) 



ire 



1 - z 



This is the same loop integral and discontinuity as in figure 20(a). For the crossed graph with 
the external lines interchanged, we can again change ^ — >• — £ to obtain the same expression 
up to d{£~— uj)6{i~^+ 6+) — )■ 6{£~ + uj)6{i~^ — b^), which leads to a vanishing discontinuity as 
before. Expanding in e, we obtain 



g/q ^ ' ^ 27r 



+m 



1 



t 



2^0 \ ,,2 1^ 91 



+ ln- 



1 - z 



+ e{i - z)z 



(A.35) 



with Pgq{z) as in eq. (A. 16). The matrix element with external antiquarks gives the same 
result, so the mixing contributions from quarks and antiquarks are identical. 

The renormalized matrix elements Bg/j(t, z, ^) are given in terms of the bare ones as [68] 



B 



bare 
9/3 



(t, z) = j dt' Z%{t - t', /i) Bg/j{t', z, /i) 



(A.36) 



Since j^{t,fi) = —fi-^Z^{t,fj,), the NLO counterterm follows from the one-loop anomalous 
dimension in eq. (A. 5) with rcusp(as) = CaolsIt^ and 7^(as) = /3oas/(27r) (see appendix B), 



Z|(t,^) = 5(t) + 



a.s(/^) 
27r 



ICa 



1 . \\ ^ ( t 
-^b[t) 2'^o( — 



+ ^/?o5(t)!>. 



(A.37) 



Expanding eq. (A.36) to NLO and using the tree-level result for Bg/j{t, z, /j,) in eq. (A. 11), 
we find 



B 



bare(l) / 



'{t, z) = Zf\t, ^^) 6g, 5(1 -z) + B^it, z,fi). (A.38) 
Note that the mixing contribution Bg/^ is UV finite and not renormalized. From eq. (A.37) 
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and our results in eqs. (A. 32) and (A. 35) we have 



27r 



2Ca 



1 



27r 



^(z)<!--<5(t)P3,(z) 



Pggiz) + 5(t) [Pg,(z) lu + ^(l - z)z 



2 9'1\ 



(A.39) 



where ah 1/e divergences are now IR divergences. Subtracting the one-loop PDF matrix 
elements in eq. (A. 18) according to eq. (A. 13), we see that the IR divergences precisely cancel 
between Bgjj and fg/j, and the finite terms in eq. (A.39) determine the one-loop matching 
coefficients Igj{t,z,fi) in eq. (2.14). 

To compare our expression for Ig^ {t, z, fi) with the result of ref. [69] , we take the moments 

Cff (z,iV) = j\yy'^l(]^[M\l-y),z,^] = dt (l - ^)\(i)(t, ^) 



2lT 



e{z) 



In^ 



(^),:-,-.(^)p„W 



^ , ,2(1 -z + z"^)^ 
+ £i(l - z)^ - Pgg{z) In z 



+ 



1 

nJ ' 



(A.40) 



where we used the approximation for the harmonic numbers at large 



N 
i=l 



(A.41) 



Our result in eq. (A.40) agrees with eq. (68) of ref. [69] except for a constant term — a^Cyi/vr (5(1- 
z)^V8- 

A. 4 Integrals, Discontinuities, and Plus Distributions 

Here we list various loop integrals that are needed for the calculation in appendix A. 3: 



1 



(27r)'^-2 (^-5 ^ 
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(A.42) 



-48- 



To calculate the discontinuities of the various graphs we need the relations 

sinvre 9{t) 
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The plus distributions are defined as 
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sinvre 9(t) 



Cn{x) 



' ' 7re(l + e) 
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(A.44) 



satisfying the boundary condition dz Cn{z) = 0. We also need the distribution identity, 
1 1 



{l-z) 



l+e 



-5{l -z) + £o(l -z)- eCi{l -z) + 0{e') . 



(A.45) 



B Perturbative Results 



B.l Hard Function 



The functions F^'^^ (z) and F^^^ (z) which encode the mt dependence of the hard Wilson coef- 
ficient in eq. (2.7) are given by 



_ 3 3 
~ 2z~ 2z 



1 

1 - - 

z 



\ arcsin^(-^) , 



< 2 < 1, 



Jn2[-i(Vi+^/^^)], z>l, 
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— / ^ 2; ^ z z I \J A 

\ 45 4725 1134 65488500 / 

/ 307 25813 2 3055907 . 659504801 ^ 

+ -3 + z + z^ + z^ + / Cf + Oi: 

V 90 18900 3969000 1309770000 / ^ 



z^). (B.l) 



Here, F^^) (z) gives the well-known mt dependence of the leading-order gg ^ H cross section 
given by the virtual top-quark loop. The full analytic mt dependence of the virtual two-loop 
corrections to gg ^ H in terms of harmonic polylogarithms were obtained in refs. [76, 77]. 
Since the corresponding exact expression for FW(z) is very long, we use the results expanded 
in mfj/m^ from ref. [79]. The additional mt dependence coming from F^^\z) is small and 
the expansion is converging very quickly, so the expanded result is completely sufficient for 
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practical purposes. For completeness we also give the leading terms in F^'^\z), 



F^'^\z) = (7Cl + IICaCf - GCfPo) ln(-4z - iO) + 



419 Tvr^ vr^ 



44C3 Cl 



+ 



217 



2 +"&j^-^'^''+V108 



/2255 5tt^ 



+ 



12 



+ 



23C3 



CaI^o — -CaTf 





27^2 



41 

y 
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(B.2) 



The first few higher-order terms in z can be extracted from the results in refs. [78, 79]. 

Alternatively to the one-step matching we use, one can perform a two-step matching by 
treating mu ^ w-t as done in refs. [49, 51, 59, 85]. In this case, one first integrates out the 
top quark at the scale — "^t and then matches from QCD onto SCET at ^ih — rnn- In 
this way one neglects power corrections of 0{mH /mt) but in turn the running between fit 
and fiH allows one to sum logarithms of mH/mt- We first integrate out the top loop to get 
the effective Hamiltonian [12, 13] 



K 



eS 
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1 + ^^{5Ca-^Cf) 
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(B.3) 



where as before as{fit) is evaluated for Uf = 5 flavors and at higher order Ci^mt, fit) contains 
logarithms of nit/ fit- In the second step, we integrate out hard off-shell modes by matching 
onto the SCET operator in eq. (2.5), 



1 + 



o^sifJ-H) ^ / , 9 -g^ - iO 
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The SCET matching coefficient C2{q'^ , fJ-n) corresponds to the IR finite parts of the gluon 
form factor. Combining the two matching steps we obtain 



CggH{mt,q^ ,fi) = Ci{mt, fi) C2iq^ , 



as{n)< 1 + 



Qs(/i) 
47r 



In^— ^-^ +5+ — Ca 



3Cf 



(B.5) 



In the second line we expanded both coefficients at the same scale fi, which reproduces the 
rrif — 7- oo limit of eq. (2.7). At NNLO, this gives the leading terms of F^'^^ (z) in eq. (B.2). The 
strict nit — )• oo limit is formally necessary for a consistent two step matching. Note that in the 
nit —7- oo limit that is often used in the literature, the full mt dependence of the leading-order 
cross section in the overall factor F^^\z) in eq. (2.7) is included, which then turns out to be 
a very good approximation for the total cross section [14, 18, 19, 97]. Considering the virtual 
corrections only, this way of defining the — )• oo limit is equivalent to the one-step matching 
in eq. (2.7) where one neglects the weak mt dependence of F^^\z) and using F^^\o). 
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B.2 Beam Function 

(2) 

We obtain the ^B-dependent terms in the two-loop contribution in eq. (2.13) by pertur- 
batively solving the two- loop RGE of the matching coefficients Igj, given by [68] 



dz' 
z' 



Xg,(t, z,fi)=J2 I dt'^^gk{t - t', ^, /i) [7^(t', Skj6{l - z') - 6{t')jl^{z', fi) 
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where the 7^(t, ;u) is the anomalous dimensions of the beam function in eq. (A. 5) and Tkji^^ f^) 
that of the PDFs. We obtain 

I^fit,z,i^B) = ^£3(^)8C7l<5(l -z) + -^C2(^) [UClPyyiz) - 2C^/3o5(l - z)] 



f^B ^f^B 



t 



LV3 



vr 



5{l~z)+2Zg^'){z) + {P,g®Pgg){z) 



y5(l -Z)- Pggiz)] + SCFTprifiPgg ^ Pgg){z] 



H T^oi—T 



4C^ 



7 . vr^ 



9 



- + 8C3U(1 -Z)- -Pgg{z) + (X(l/) 



+ CA/3o[(-f + y)5(l-^)-2XS'^)(. 



+ CFTFnf[A5{l -z) + 8(X(i'^) ® P,,)(z)] + 4Pg«(^)| + 45(0 ^f^^) , 
xg)(t,z,^B) = ^£2(^)l2CACFP5,(z) + ^A(^)|4C|.(P<;g^5P,g)(z) 

/^B /^B ^^B f^B I 



f^B 



-^P,,(z) + (X(i'^)0P,,)(z) 



(B.7) 



+ 4C|(l£'^) ® P,,)(z) - iCF(3o4Y\z) + AP^f{z) + 4<5(t)x(2.5)(^) 



(0 A"! 

The /iB-independent terms X (z) require the two-loop calculation of gluon beam function 



which is not yet available. The one-loop functions x'^g-'^\z) entering in the above expressions 



were already given in eq. (2.15). The quark and gluon splitting functions are 
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The two- loop splitting functions were calculated in refs. [98, 99] and are given by [99] 
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The convolution of two splitting functions is defined as (the index j is not summed over) 
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and analogously for the convolution {llj'^^ Pjk)- The necessary convolutions are 
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2 - - + (2 - z) In z + ^ ln(l - z) , 



-31 + 24z + 3z2 +4^3 4(l + z + z2) 2(2 - 2z + ^2) 

Inz -\ ln(l — z) 



3z 



27r2 



(Pgg Pgg){z) = 8£l(l " Z) " — <5 ( 1 " z) + 



z z 
4(-ll + 9z - 9z2 + 11^3) 



4(-l-3z2+ 423-^4) 

+ -, ^ In z + 



3z 

Kl-2z + z2-z3) 



{X^^/^ ® Pgg){z) 



(X(i'^)®P,,)(z) 



z(l — z) z 
13 + 12z + 6z2 - 5z3 o 1 + z -4 + 9z2 + 4.z^ 



9z 



+ 7r^ 



+ 



3z 



ln(l - z) , 
In z 



4 + 3z - 3z2 - 4z3 , , , , , n 
+ ln(l - z) - (1 + z) ln2 z - 2(1 + z)Li2(z) , 

oZ 

5-4z + 2z2 4 + 6z-3z2 -2 - z , 4 + 3z , , 

+ vr^ \ — In z H — ln(l - z) 



2z 



6z 
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+ z + ^ ln(l - z) [ln(l - z) - In z] + (z - 2)^2(2) , 

.^ „ , , , 21 - 26z + 2 - 6z - 9 - 30z - 9^^ - 42^ 

(^S'*' « = + Si + 1- 

^-31 + 24,- + 3.^ + 4.3 2(1 + ,- + .^)^^,^ 

2 — 2z -I- 

+ -— ln(l - z)nn(l - z) - 2 In zl + (8 + 2z)Li2(z) , 

z 

(4^''^ ® P,,)i^) = 6^:2(1 - - T^'Coil -z) + 4C35(1 - z) + + 

n -3 + 2z - 7z2 + 3z3 2(11 - 21z + Gz^ - 22z3) ^ 

+ TT^ h ^ ^ In z 

3z 3z 

4(-ll + 9z - 9z2 + llz^) ^ ^ ^ 2(1 + 3z2 - 4z3 + z"^) , n 

+ - 5; -'""-''+ — ^'"^ 

z(l - z) 

+ 6(1-2^ + ^'-^') i^^2(-^ _ ^) ^ ^ ^)Li2(z) . (B.ll) 
z 

B.3 Renormalization Group Evolution 

The hard Wilson coefficient satisfies the RGE 

7^(9',/^) = r?usp[«.(/^)] In^^^^ + 7|fK(/x)] , (B.12) 

A* 

whose solution gives eq. (2.10). The beam function RGE is given in eq. (A. 4). The soft 
function satisfies an analogous RGE, 



Ai^5'f(A;,^) = j dk' -fl{k - k' , fi) S'^jf{k',n) , 

7|(A;,/i) =4r?,3p[a,(/i)]^/:o(^) + 7|[as(/i)] <5(A:) , (B.13) 

whose solution yields eq. (2.24). The functions i^p(^0;A*); ^r(/^0;/^)) K^{^q,^) required for 
the RGE solutions of the hard, beam, and soft functions in section 2 are defined as 



Kh{l^o,^^)= ^r^r* (a,) / T^TTT' ^r(/^o,/i) = / ^rA-T* (a^). 



i^,(A^o,/u)= / ^^^^7(a,). (B.14) 
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Expanding the beta function and anomalous dimensions in powers of ag, 



/3(a,) = -2a, 



n=0 



47r/ 



'V47r 



7(0.) = XI' 



47r 



5 -'^cusp('^s) ~ ^ ^ -^n ( 

re=0 

(B.15) 

their expHcit expressions at NNLL are (suppressing the superscript i on K^, r/p and F^), 

^r(/xo,/.) = -^|44(l---lnr) + (^-|)(l-r + lnO + ^ln^r 
4/?o l"s(/Uo) V r / \To /3o/ 2/3o 



+ 



0) 



Air 
£2 

ro (^0^0 



2 



1 



+ In r I + 



/3oro 1^1 



(1 — r + r Inr) 



/3iri\ (1-0' 



^r(/io,/^) 



'2/3o 



Inr + 



To /3o 



a 



70 
2/3o 



+ 

Inr + 



47r 



(r-l) 



/3iri , /32\r2-l 



+ 



(r-l) 



(B.16) 



167r2 /3oro /3o' /^o, 

4vr V7o /5o, 

Here, r = as(/i)/as(/Uo) and the running coupling at the scale fj. is given in terms of that at 
the reference scale /io by the three-loop expression 



X 



167r2 



/32 



/3o 



1\ (3ff\nX 1 



/32 V X 



(B.17) 



where X = 1 + as(/xo)/3o ln(/i//io)/(2vr). At NNLL, we require the three-loop (3 function 
and gluon cusp anomalous dimension and the two- loop non-cusp anomalous dimensions. The 
coefficients of the MS /3 function to three loops are [100, 101] 

Po = y C'a - - 7> n/ , 

/3i = yCl- (yC-A +4CF)TFnf, 



f32 



2857 
~5A 



Ca+ [Cp TTT CfCa 



1415 



11 



79 



CI ) 2ri. + ( — CiT + — CaJ nj . (B.18) 



18 ""^ 54 

The coefficients of the gluon cusp anomalous dimension in MS are [102] 

rg = 4CA, 



rf = AC A 



V 6 



27 ^ 45 ^ 3 



27 27 



Ca Tf uf 



55 



16 



+ (-^ + 16C3)CFTi.n/--^r|n^ 



27 



(B.19) 
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and are equal to C a/ Cf^q 12- The anomalous dimension coefficients for the hard Wilson 
coefficient to three loops are [59, 103, 104] 

7^o = -2/5o, 

7^i = (-^ + 4C3)cl+(-f + ^)c7./3o-2/3,, 
g _( 60875 6347r2 Svr^ 1972C3 4071^(3 ^3 

The non-cusp anomalous dimension for the gluon beam function is equal to that of the gluon 
jet function which is given by 7j(as) = — 27|^(as) — 7j(as). Here, 7j(as) is the coefficient 
of the 5{1 — z) in the gluon PDF anomalous dimension which is known to three loops [102]. 
The resulting coefficients to three loops are 

71.0 = 2/3o, 

7|, = (^-32a)ci+(^-^)c.A + 2A, 

,49373 _ 944^ _ 16^ _ 45206 1^ 

'^2 V 81 81 5 9 9 / ^ 

61 73 3767r2 ^ 13^ ^ 280(3 ^ 986 lOvr^ _ 56(3 ^ „2 



27 81 5 9 / ^' V 81 

2 0^4 



(B.21) 



/1765 27T' 87r4 304(3 ^ 
+ - — - 1^ - ^r^^^ + ^^^- 

At NNLL we only need 'jfj and 7^ at two loops. The three-loop coefficients are given for 
completeness. The result in eq. (B.21) agrees with that given in ref. [105] for the gluon jet 
function.^ The non-cusp anomalous dimension of the gluon beam-thrust soft function is given 
by 7s("^) = -27|f(as) - 2jf^{o:s)- 

B.4 Singular Fixed-Order NLO and NNLO Coefficients 

In this Appendix we give our results for the perturbative coefficients Cij entering in the fixed- 
order cross section in eq. (2.30). The coefficients for quarks and antiquarks are related as 
follows: 

Cgcj{Za, Zli, r, Y, /i) = Cgq{Za, Zb, T, Y, /i) , Cqg{Za, Z^, T, Y, fl) = Cqg{Za, Zf,, T, Y, /i) , 
Cqg{Za, Zf,, T, Y, /i) = Cgq{zi), Za, T, —Y, /i) , 

Cqq'iZa, ^fe, T, Y, fl) = C|g,(Za, Zb, T, Y, fl) = C~,{Za, Zb, T, Y, fi) = Cgg,{Za, Zb, T, Y, n) . (B.22) 



* Apart from a typo in ref. [105] where one of the terms in the C^t^/ contribution is missing a vr^ 
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The first line follows from charge-conjugation invariance and the second line from parity 
invariance of QCD. The last line is only true because we limit ourselves to contributions 
from gluon beam functions. Consequently, we only need to consider Cgg, Cgg, and Cgq'. The 
small contributions involving quark beam functions that we neglect correspond to qq ^ H 
production, where the qq pair couples to the Higgs either directly or indirectly through a 
two-loop qq ^ gg ^ H diagram. 

As in eq. (2.31) we split the coefficients into singular and nonsingular parts, where the 
singular coefficients are further decomposed as in eq. (2.32), 

C^Za, Zfe, T, y, ^l) = Cr^{Za, Zb, Y, fl) 6{t) + ^ C^jiZa, Zfo, Y, fl) Ck{T) . (B.23) 

fc>0 

Each of the coefficients has an expansion in Og, which we write as {x = s or x = k > —1) 

corresponding to the LO, NLO, and NNLO contributions. (The overall of the Higgs cross 
section is factored out in eq. (2.30).) At leading order, the only nonzero coefficient is 

Cf^\za,z^„T,Y,^Ji)=5{T)5{l-Za)6{l-Zb). (B.25) 

The singular NLO terms are fully contained in the resummed result at NNLL. Hence, the 
simplest way to obtain them is to set fin = fJ'B = l-i-s = fJ' eq. (2.2), which eliminates the 
evolution factors, and expand it in as(/x). Using the NLO results for the hard, beam, and 
soft functions, we find 

Cf^\za,Zb,T,Y,fx) = CASil-Za)e{zb)\-2Ci{T)5{l-Zb) + Co{T)Pggizb) 



/27r2 . fW 
I + Y^ + 



6(1 - z,) - (y + 2L^)Pgg{z,) +Xg^'\z,) I 



V 3 2Ca' 

+ {Za O Zfe,y ^ -Y) , 
Cfpiza, Zb, r, Y, fx) = Cf 6{1 - Za) e{zb){Co(.r)Pgg{zb) 

+ 5(r) [4Y\z,) - (y + 2L^)Pgg{z,)] } . (B.26) 

Here = ln{fi/m}{) and F^^^ = F^^)[m|^/(4mj)]. The terms in the NLO coefficients 
precisely cancel the fi dependence in the leading-order cross section coming from eq. (B.25). 

k(2) 

The NNLO coefficients C^j for A; > are reproduced by the resummed result at NNLL. 
Once we include the //-dependent NNLO terms in the hard, beam, and soft functions (which 
are determined from the NNLL resummation), we can obtain the c'^j'^^ as before by setting 

I^H = Mb = /^s = M ill eq. (2.2) and expanding to 0{a1). For C^j^^"^^ we are only able to 
obtain some parts analytically. We write it as [see eq. (2.38)] 

Cr^^'^^Hza, Zb, y, fl) = cJAZa, Z^, Y) + cf,(Za, Zb, Y, fi) + cjf (z„ z^, Y) . (B.27) 
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The first contribution denotes the vr^ terms, which we get by taking fin = M = nT-H in 
the resummed result (as opposed to taking = —irrifj, which resums them into the hard 
evolution factor). The second contribution contains the /x-dependent terms proportional to 
= ln(/i/m/f), which we are able to obtain by requiring that they cancel the fj, depen- 
dence of the singular NLO result. Obtaining an analytic expression for the remaining piece, 
(f-^^{za, Zb,Y), requires the complete two-loop hard, beam, and soft functions. Its contribu- 
tion to the cross section after convolution with the PDFs and integrated over Y is extracted 
numerically from the fixed-order NNLO cross section as discussed in section 2.4. 

s(2) 

For the nonzero coefficients of Cgg we find 

C^g^^\za, Zb, Y, ll) = 8Cl J(l - Za)6{l - Zb) , 

6-2(2) (z„Zb,y,/x) ^ _ Za)[-QClPgg{zb) + ^Ca/3o <5(1 -Zb)} + {Za ^ Zb) , 

Clf\za, Zb,Y,fi) = 6(1 - - yvr^ - AY^)6{1 - Zb) + 4(y + 2L^)Pgg{zb) 

+ {Pgg ® P99)i^b) - ^'i'/Kzb)] - 2CaF(i)5(1 - Zb) 

- CApa[{^ + 2L^)<5(1 - Zb) + \Pgg{Zb)] + 2CFTFnf{Pgg Pgg){Zb)] 
+ ClPgg{Za)Pgg{zb) + {Za ^ Zb,Y ^ -Y) , 
C^gP{Za,Zb,Y,fi) = 6{1 - Za){cl[{l + 5C3)<5(1 - Zb) + 2(7r2 + Y^)Pgg{zb) 

- (y + 2L,){Pgg Pgg){Zb) + Pgg){Zb) 

.2 



^Ca^o [(2 + ^ + y2)<^(l - Zb) -{Y + 2L^)Pgg{zb) +lg''\zb) 



+ 2CFTFnf[^6il - Zb) -{Y + 2L^){Pgg ® Pqg){zb) + (X^J'^) ® Pqg){zh) 

Tg'^Kza) - 2L^Pgg{Za)\ PggiZb) + {Za ^ Zb,Y ^ -Y) . (B.28) 



The two- loop splitting function, Pgl\z) is given in eq. (B.9). The convolutions {Pij ® Pjk){z) 
and (Xj-j^'"^^ (g) PjjS)[z) are defined in eq. (B.IO) and their analytical expressions are given in 
eq. (B.ll). For the //-dependent and vr^ terms of Cg^^"^^ we get 



cUZa,Zb,Y,f,) = -2L^5{1 - Za)\ C\ 



^2Y'']Pgg{zb 



(y + L^){Pgg ® Pgg){Zb) + 1^^^ ^ Pgg{Zb) 



(B.29) 



+ Y')6{1 - Zb) - (y + L,)Pgg{Zb) 



V 3 



-gg 



[Zb 



+ 2CFTFnf -(y + L^){Pgg Pgg)iZb) +I^gY^ ® Pggi 



Zb 



i(/3i + f(^)/3o)5(1 - Zb) + CAF^'^Pggizb) + P^l\zb)} 
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2Cl - L^PggiZa)] Pggizb) + {Za ^ Zb,Y ^ -Y) , 



99 



{Za,Z,,Y) = T,H{1 - Za){cl[[- + ^ +Y^)6{1 - Z,) - Y Pgg{zb) + lg''\z,) 

+ ^C^/3o <5(1 - Zb) + -Zb)} + {za o zfe, y ^ -Y) . 



For Cgq^^ we find 



Clf\za,Zk,Y,^i)=Q, 



(B.30) 



Cf^\za,Zi„Y,^i) 



-QCACF6{l-Za)Pgq{zi,), 

Clf\za, Zfe, y, ^i) = 5(1 - Za)[CACF [4(y + 2L^)Pgg{zb) + {Pgg ® Pg,){zb) - 4Z(J'^)(zb) 

+ CUPgq ® Pgg)iZb) " /^O ( ^6 ) } + 2C aC p PggiZa) PggiZb) , 

Clf\Za, Zfe, y, = - Za){CACF [2{^^ + Y^)Pgg{Zb) - {Y + 2L^){Pgg ® Pgg){Zb) 

+ {I^l''^ Pgg){Zb) 

+ C| [-{Y + 2L^){Pg, ® P,,){Zb) + (X(J'^) ® P,g){Zb) 

+ CFfio[{Y + 2L^)Pgg{zb)-I^gy\zb)\ +CFF^'^Pgg{zb) + 
+ CaCf 



Pgg{Za)lf/\zb)+X^/\za)Pg,{Zb)-AL^Pgg{Za)Pg,{zb) 



-99 



The two-loop splitting function Pgq {z) is given in eq. (B.9). For Cgq^'^' , the /i-dependent 
and TT^ terms are 



1(2) 



(Za,2b,y,/i) = -2L^(5(1 



){CaCf 



(— + )p,,(z,) - (y + L^)(P33 ® p,,)(z6) 



+ (Xg''5) ® Pg,)(zfc)] + [-(y + L^)(P,, ^ P,,)(z,) + (X(i'^) ® P,,)(z,) 
+ CFf^o [{Y + L^)Pg,{zb) - X^^/\zb)\ + CFF^^^Pg^izb) + P^']^{zb)} 

- 2CACFL^Pgg{Za)X'^^/\zb) +lg'^\za)Pgq{Zb) " 2L ^ Pgg{Za)Pgq{Zb) 



c^g^{za,Zb,Y) = tt^CaCf5{1 - Za)\-YPgq{zb) +X'i'/\zb) 



(B.31) 



Finally, for cf^^ we find 

.3(2), 



C;';' [Za, Zb, Y, fl) = C^^^ {Za, Zb, Y, f,) = , 

Clgf^ (Za,Zb, y, fl) = 2Cp Pgg{Za)Pgq' (zb) , 
-,0(2) 



Clf{Za,Zb,Y,^i) = Cl\Pg,{Za)X^'^P{Zb) + X^^/\za)Pgq' {Zb) - ^ Pg,{Za)Pg,> {Zb) 



(B.32) 
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There is no contribution to C^^/^^^ from vr^ summation, c^^, = 0. The /i-dependent terms are 



qq' 



C^.{Za,Zb,Y,fl) = -2ClL^ Pgq{Za)lllf\zt)+I^YHza)Pgq'{Zb) " 2L^ Pg,{Za)Pgq' (zb) 



11 



91 



(B.33) 
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